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Zusammenfassung 

Die vorlicgcndc Arbeit bcschaftigt sicli iriit (k;r supcrsyinnic'trischcn Vc^rallgc- 
meinerung von Matrix- und Eigenwertmodellen. Nach einer kurzen Einfiihrung in 
das Hermitesche Ein-Matrix-Modell wenden wir uns dem c = — 2 Matrix-Modell 
zu. In seiner FormTilieriing durch cin matrixwcrtiges Snperfcld ist dieses Modell 
invariant unter Supersymmetrietransformationen auf Matrixebene. Wir zeigen die 
Existenz einer Nicolai-Abbildung dieses Modells auf ein freies Hermitesches Matrix- 
Modell und diskuticrcn seine diagrammatische Entwicklung. Korrclationsfunktio- 
nen fiir quartische Potentiale und beliebigen Genus werden berechnet, welche die 
Stringsuszeptibilitat von c = — 2 Liouville-Theorie im Skalenlimes aufweisen. Wir 
zeigen auf, wie sich diese Ergebnisse zum Zahlen supersymmetrischer Graphen ver- 
wenden lassen. 

Daraufhin studieren wir das Supereigenwertmodell, der bis heute einzige erfolgre- 
iche diskretc Zugang zur Quantisierung von 2d Supergravitation. Das Modell wird 
in einer superkonformen Feldtheorie Formulierung durch Forderung von Super- 
Virasoro Bedingungen konstruiert. Die vollstandige Losung wird mit Hilfe der Mo- 
nientcnmcthode hcrgeleitet. Dicsc ermoglicht die Bcrcchnung der frcicn Energie 
und aller Multi-Loop-Korrelatoren auf beliebigem Genus und fiir allgemeine Po- 
tentiale. Die Losung wird im diskreten Fall und im Doppelskalenlimes prasentiert. 
Explizite Rcsultate bis Genus zwei werden angegeben. 

Es folgt eine Diskussion der supersymmetrischen Verallgemeinerung des externen 
Feld Problems. Die diskreten Super-Miwa-Transformationen des Supereigenwert- 
modells konnen im Eigenwert und im Matrixfall angegeben werden. Eigenschaften 
von externen Supereigenwertmodellen werden diskutiert, die genaue Form des zum 
gewohnlichen Modell korrespondierenden externen Supereigenwertmodells konnte 
jedoch noch nicht hergeleitet werden. 
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SUPERSYMMETRIC GENERALIZATIONS OF MATRIX MODELS 



Abstract 

In this thesis generalizations of matrix and eigenvalue models involving supersym- 
metry are discussed. Following a brief review of the Hermitian one matrix model, 
the c = — 2 matrix model is considered. Built from a matrix valued superfield this 
model displays supersymmetry on the matrix level. We stress the emergence of a 
Nicolai-map of this model to a free Hermitian matrix model and study its diagram- 
matic expansion in detail. Correlation functions for quartic potentials on arbitrary 
genus are computed, reproducing the string susceptibility of c = —2 Liouville theory 
in the scaling limit. The results may be used to perform a counting of supersym- 
metric graphs. 

We then turn to the supereigenvalue model, today's only successful discrete ap- 
proach to 2d quantum supergravity. The model is constructed in a superconformal 
field theory formulation by imposing the super-Virasoro constraints. The complete 
solution of the model is given in the moment description, allowing the calculation 
of the free energy and the multi-loop correlators on arbitrary genus and for general 
potentials. The solution is presented in the discrete case and in the double scaling 
limit. Explicit results up to genus two are stated. 

Finally the supersymmetric generalization of the external field problem is addressed. 
We state the discrete super-Miwa transformations of the supereigenvalue model on 
the eigenvalue and matrix level. Properties of external supereigenvalue models are 
discussed, although the model corresponding to the ordinary supereigenvalue model 
could not be identified so far. 
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Introduction 



Random matrices have been studied in physics since the work of E. Wigner in 
the 1950's. Initially proposed as an effective model for higher excitations in nuclei, 
they have found numerous applications in various fields throughout the years. In 
high-energy and mathematical physics matrix models experienced a great renais- 
sance following the discovery of their relevance for the quantization of two dimen- 
sional gravity and bosonic string theory in 1985. Only recently the supersymmetric 
generalization of matrix models has been addressed. This thesis reports on the con- 
struction, solution and interpretation of such generalized models, focusing on the 
work of the author. 

Matrix models have the appealing property that they are exactly solvable in the 
limit of infinite matrix size N. Subleading corrections to arbitrary order in 1/iV may 
be computed by iterative means, while results for finite N are few. The analysis of 
the diagrammatic expansion of the Hermitian matrix model reveals that solving the 
model is equivalent to performing the sum over random equilateral triangulations 
of two-dimensional surfaces of fixed genus. 

String theory on the other hand is today's most popular candidate for a unification 
of quantum field theory and gravity. It is built on the attractive idea of replacing 
particles by string-like one-dimensional objects, which sweep out a two-dimensional 
"world-sheet" as they evolve in time. According to Polyakov this theory is to be 
interpreted as two-dimensional quantum gravity, in which the world-sheet replaces 
two-dimensional space-time, and where the coordinates of the string form matter 
fields coupled to the two-dimensional gravity theory. Interactions correspond to 
topology changes and are encoded in the genus of the world-sheet, which may be 
viewed as an "inflated" Feynman graph. Thus the task in a path-integral quantiza- 
tion is to perform an integral over the two dimensional geometries of the world-sheet 
and sum over their genera. 

Precisely this integration may be performed with the help of the matrix model 
by taking a "continuum limit" of the above-mentioned exact solution, i.e. when the 
triangles of the discretized surfaces become dense. However, this rather unconven- 
tional method of performing a path integral only works for toy models of bosonic 
strings with world-sheets of Euclidean signature living in Z? = c < 1 dimensions. 
Extending these techniques to dimensions greater than one is a hard problem still 
unsolved. Due to the exact solvability of matrix models new non-perturbative tech- 
niques were introduced to these low-dimensional string theories. We will review the 



Hermitian matrix model and its correspondence to random surfaces in chapter I. 

Of potential phenomenological relevance, however, are superstring theories, which 
in contrast to the bosonic string have bosons and fermions in their spectrum. The 
Neveu-Schwarz-Ramond superstring is to be interpreted as two dimensional super- 
gravity coupled to a matter action of superfields. It is thus very desirable to find a 
supersymmetric generalization of the matrix model technique in order to perform 
the corresponding path integral over super-geometries. There are (at least) three 
different ways to proceed. 

In a geometrical approach one puts superstrings on a random lattice and then 
tries to establish a map to an adequate supersymmetric matrix model. But due 
to the known problems with supersymmetry and fermions in lattice field theory 
this seems to be rather hopeless. However, the alternative superstring theory of 
Green and Schwarz does not suffer from these problems, as here the fields on the 
world-sheet remain bosonic and supersymmetry is only present on the external 
space-time level. A lattice version of this superstring theory has been formulated 
and studied. Still, the relation to a solvable supersymmetric matrix model is an 
open problem (for more details see refs. [1]). Even classically the Green-Schwarz 
string is supersymmetric only in _D > 4 dimensions, a domain that seems to be 
hardly reachable for solvable matrix models. 

Alternatively one can consider supersymmetric generalizations of matrix models 
and hope that their critical behaviour reveals properties of two-dimensional super- 
gravity. However, no example is known in the literature where this is the case. 
Nevertheless, the analysis of these models is interesting in its own right from the 
viewpoint of random matrices. We shall therefore study the c = —2 matrix model 
in chapter II, a model displaying supersymmetry on the level of matrices. In its 
"continuum limit" this model describes the propagation of bosonic strings in minus 
two-dimensional space-time. 

The final and successful approach seeks for a generalization of the integrable 
structure of the Hermitian matrix model. This structure at the heart of the Her- 
mitian matrix model resides in a set of Virasoro constraints. It turns out that, by 
imposing a generalization in the form of super- Virasoro constraints, one obtains 
the supereigenvalue model of Alvarez-Gaume, Itoyama, Manes and Zadra, which 
generalizes the eigenvalue formulation of the Hermitian matrix model. Due to the 
eigenvalue (rather than matrix) character of this model there is no geometrical 
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Introduction 

interpretation in terms of triangulated super-Riemann surfaces at hand. The su- 
pereigenvalue model is exactly solvable for general polynomial potentials in the limit 
of an infinite number N of eigenvalues. Moreover, all subleading corrections in 
are determined through an iterative process. We will describe this solution in chap- 
ter III in the "discrete" and "continuum" cases. The continuum results reveal that 
the supereigenvalue model indeed describes the coupling of minimal superconformal 
theories to two dimensional quantum supergravity. 

As we shall see, the supereigenvalue model has a large number of similarities 
to the Hermitian matrix model. This motivated the author to study the external 
field problem in this context as well. In chapter IV we review the external Hermi- 
tian matrix model and present some preliminary results of a generalization to the 
supersymmetric case. 

The presentation of this thesis is intended to be self contained, only basic knowl- 
edge of conformal field theory is assumed. For a review on this topic see e.g. ref. [2]. 
There is a large number of reviews on matrix models and two-dimensional quantum 
gravity, e.g. refs. [3], but only the recent ref. [4] contains the supereigenvalue model. 

The parts of this thesis containing the solution of the supereigenvalue model have 
been published prior ly in refs. [5, 6]. 
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Chapter I: 

The Hermitian Matrix Model 



Generally speaking matrix models are quantum field theories where the field is 
a. N X N real or complex matrix M{x). We shall consider the simple case of the 
Hermitian one matrix model in Z? = dimensions, which due to its simplicity 
in exactly solvable in the limit of infinite matrix size N. This model has been 
intensively studied in the literature. In the following we give a basic introduction 
to the model, as well as a review on a collection of more detailed aspects which will 
be relevant for the generalization to the supersynimetric case. 



1. The Model 

The Hermitian one matrix model is defined by the partition function 



Zn [gk] 



VM exp 



-NTiV{M) 



(1.1.1) 



where Mij is a N x N Hermitian matrix. F[gk] is the free energy. The measure for 
the "path integral" of this zero dimensional theory is given by 



(1.1.2) 



and we consider the most general polynomial matrix potential with coupling con- 
stants gk 



(1.1.3) 



fc=0 



Note the U{N) invariance M^WMU of the "action" TrV{M). This 
may be used to diagonalize the Hermitian matrix M ^ D U , where 
D = diag(Ai, . . . , Aat) is the diagonal matrix. By performing the change of vari- 
ables from M to Xi and U, the integral over the unitary group factors out and we 
are left with the integration over the eigenvalues Xi . The Jacobian of this transfor- 
mation is YliKji^i - ^if- 



This may be seen by considering the norm of the infinitesimal variation of M = DU 
I 5M\'^ = ^ SMij SMji = Tr {5Mf 



■.Tt(^-UUuu'' DU + 5DU + U'' D5U 



= Tr {SDf - 2i Tr [5D, D]5u + 2 Tr {-Su DSuD+ {Suf ), 

where we have introduced Su = i SU = Su^ . The second term in the last expression vanishes as 
SD and D arc diagonal. We then find 

I |2 = '^{S\,f + ^(A, - \,f \Su,,\\ 

Note that the independent variables are the variation of the eigenvalues S\i and 9ieSuij , 3mSuij 
for i < j. From the last equation we obtain the Jacobian of M t~> Xi,U which is y/G where G is 
the metric tensor, explicitly 



The Hermitian one matrix model (1.1.1) may then be written in the eigenvalue 
representation 



N 



N 



ZN[9k] = CN / ([[dX,) Y[{X, - X,y exp ~N ^ V{X,) 



(1.1.4) 



i<j 



where cat is the volume of the U{N) group. Note the relation of the Jacobian A^(A) 
to the van der Monde determinant A(A) = (-)^(^-i)/2 Det (A^"^). 

Loop Insertion Operator 



Expectation values in the Hermitian one matrix model are defined in the usual way 

as 

(1.1.5) 



'O(M)) = ^ / VMO{M) exp -NTi-V{M) 
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I: The Hermitian Matrix Model 

A similar expression holds in the eigenvalue picture. It is very convenient to work 
with the one-loop correlator 

k=0 ^ 

which acts as a generating functional for the amplitudes ( Tr Af'^ ) . Similarly we may 
define the generating functional for higher point amplitudes, the n-loop correlators 

,,rr ^ ^ ( Tr M''! . . . Tr M'^") conn. ,^ ^ 

fei,...,fc„=i Pi ■■■P" 



where "conn." refers to the connected part. The last two equations may be rewritten 
as 

W{pu...,pn)^N"-^(TT ^rT...Tr ^-—) . (1.1.8) 

\ pi — M pn — M / conn. 

A useful object to define is the loop insertion operator [8] 

S _ ^ 1 d 

because we may now obtain the rt-loop correlators from the free energy F by ap- 
plying the loop insertion operators: 

Hence once the free energy for a general potential is known, all observables may be 
calculated. The same holds true for the one-loop correlator W(pi), as with 

mP....,P«) = ^...^m.r) (1.1.11) 

all multi-loop correlators follow from it. We thus see that solving for W{p) with a 
general potential really means completely solving the Hermitian one matrix model. 
We shall see in section 4 that this is most efficiently done by considering the loop 
equations of the model. 

It is known since the work of t'Hooft [9] that the free energy F and all correlators 
admit an expansion in which may be seen by looking at the perturbative 

evaluation of the matrix model (1.1.1). 



1.1. The Model 




Figure 1: The propagator and a four point vertex of the Hermitian matrix model 



Feynman Diagrams 

The propagators of the Mij fields of eq. (1.1-1) may be represented diagrammat- 
ically by double lines each one corresponding to the separate propagation of the 
matrix indices. For Hermitian matrices the lines should be oriented in opposite 
directions (cf. Figure 1). The propagator then simply is 



' Mij Mki )o ^ —Sik Sji, 



(1.1.12) 



where the subscript denotes the average taken in the free theory, i.e. g2 = 1/2 and 
gk = for k ^ 2. With the general interactions of eq. (1.1.3) there will be three- 
point vertices, four-point vertices, etc. . . , each n-point vertex contributing a factor 
of {—gn N). Moreover each loop of internal index will yield the factor N = Su. 

Consider for instance a connected vacuum diagram built out of PCpSopagators, 
L closed loops of internal index and V3 three-point vertices, V4 four-point Yefj;ices, 
etc. . . , and let 

V = V3 + V4 + V5 + .... 



(LI. 13) 

Each loop of internal index may be considered as a face of a polyhedron, and the 
Euler relation then gives 

V-P + L = 2-2g, (1.1.14) 

in which g is the genus of the surfaces on which the polyhedron (or Feynman di- 
agram) is drawn (0 for a sphere, 1 for a torus, etc. . . ). The contribution of this 
diagram is 

N'^-^+^i-g^fH-g^f^-- - =N'-'^l[{-gkr''. (1.1.15) 



I: The Hermitian Matrix Model 

As the free energy F is the sum of aU connected diagrams we find the genus expan- 
sion 



(l.f.l6) 



where Fg is the sum of aU diagrams that can be drawn on a surface of genus g. 
The Hermitian one-matrix model can be solved rather easily in the "planar" limit 
of infinite matrix size iV ^ cxd as we shall see in section 4. Here only the planar 
diagrams of Fq survive, the Hermitian matrix model thus represents a powerful tool 
for planar graph counting [10]. 

2. Discretized Surfaces and 2d Quantum Gravity 

We learned in the previous section that the closed Feynman diagrams of the Her- 
mitian matrix model may be viewed as polyhedrons. Moreover the matrix size N 
appeared as the parameter controlling the genus of such a polyhedron. The nontriv- 
ial combinatorial problem of how many inequivalent polyhedrons of genus g with V3 
3-point vertices, V4 4-point vertices, etc. . . . exist, may be answered by expanding 
^s[5fe] of eq. (1.1.16) in the gkS. But as a polyhedron may be interpreted as a 
discrete approximation to a smooth two dimensional surface, this result may be ex- 
ploited to study an entirely different problem: The quantization of two ( ifam ensional 
Euclidean gravity! 

In two dimensions the Einstein-Hilbert action forms a topological invariant, the 
Euler characteristic of the underlying manifold. So for fixed topologies only the 
cosmological term will be dynamic. If we consider a specific manifold the action of 
2d Euclidean gravity is 

S^^x j dx^^g j dx^^gR=^xA-'?-^ (1.2.1) 

where G denotes the gravitational constant, /i the cosmological constant and A the 
area of the surface. In a path integral quantization of this theory the integration 
over the metric may be split up into separate integrals over topologies and aj; 



J VoKDiff) ^ Jo Js,,A Vol(D«if) 

where we have formally divided out the volume of the group of diffeomorphisms. In 
general the volume of the moduli space Vol(I]g.yi) — ^ 'Dg^v /voi{Diff) is difficult 
to calculate. The most progress in quantizing 2d gravity in the continuum has been 
made via the Liouville approach [11]. If we discretize the surface, on the other hand. 
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1.2. Discretized Surfaces and 2d Quantum Gravity 



it turns out that (1.2.2) is much easier to calculate, even before removing the finite 
cutoff. We consider in particular a "random triangulation" , in which the surface is 
constructed from equilateral triangles. ^ Assign the tmit area e to each triangle, 
so that the total area of this simplicial manifold built _o it of F triangles is given 
by e F. The path integral over the metric J Vgf^^ is now replaced by a sum over 
triangulations, so that the discretized partition function Zqq of eq. (1.2.2) reads 



--QG 



EE 

9=0 _F=1 



(2-2s)/Gg-peiL^ 



(1.2.3) 



where T^g^F denotes the total number of inequivalent triangulations of genus g built 
out of F triangles. 

In the above, triangles do not play an essential role and may be replaced by any 
set of polygons leading to a similar expression as in eq. (1.2.3). 

As we saw in the previous section J^g^F is calculated by the free energy of the 
Hermitian matrix model. In particular to establish contact to random triangulations 
consider the matrix model of eq. (1.1.1) with a cubic interaction, i.e. 52 = 1/2, 
gs — X and g^ ~ for fc > 3. The dual to a Feynman diagram of this model (in 
which each face, edge and vertex is associated respectively to a dual vertex, edge and 
face) is identical to a random triangulation of some orientable Riemann surface (cf. 
Figure 2). Hence there is a one-to-one correspondence between Feynman diagram 
vertices and the faces of the random triangulation. By looking at the free energy of 
this cubic matrix model 



iV^ F[X] 



g=0 V=0 



N^-^^-XfFgy, 



(1.2.4) 



we see that by formally identifying 

e^/G = N and 



-Me 



-A 
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2.5) 



in eqs. (1.2.3) and (1.2.4) the free energy of the matrix model is actually the partition 
function of discrete 2d gravity 



iV2F[.gfe]=Z^'G 



(1.2.6) 



^This constitutes the basic difference to Regge calculus, where the link lengths are the geomet- 
rical degrees of freedom. 



I: The Hermitian Matrix Model 




Figure 2: A piece of a random triangulation of a surface. Each of the triangular 
faces is dual to a three point vertex of a Hermitian matrix model. 

Needless to say that this argument goes through for general matrix potentials 
of eq. (1.1.3) as well, then describing "random polygonulations" of surfaces. In 
the continuum limit the result should not depend of the particular discretization 
used. However the different possibilities for generating vertices constitute additional 
degrees of freedom that can be realized as the coupling of 2d gravity to different 
varieties of matter. We shall see how this works in chapter HI in the context of the 
m'th multicritical points [12]. 

Continuum Limit 

In order to discuss the continuum limit we have to anticipate a result which will 
be derived in later chapters. The genus g contributions to the free energy of the 
matrix model are analytic functions of the couplings gk around the origin. Singular 
behaviour arises at a critical surface in the space of couplings. Specializing again 
to the cubic model of eq. (1.2.4) this means that by moving the constant A to its 
critical value Ac the free energy becomes singular. In the vicinity of the critical 
point one finds the scaling behaviour 

Fg[\]^iX-\,)^^-''°^^'-s), (1.2.7) 

in the case of the cubic potential the critical exponent 70 takes the value 70 = —1/2. 
We can extract the continuum limit of the random triangulation by tuning A ^ Ac. 



1.2. Discretized B+trfaces and 2d-Ouantum Gravity 



This is because the atftm 



age value of the area of the surface is g[vep by 
e{V) 



d 



eX-luF, 



A- Ac 



(1.2.. 



er of vertices V 



appearmg as 
ith respect to 



(recall that the area is proportional to the nu 
the power of A associated to each graph which has to be averagedU-w 
Fg [A] , as the free energy of the matrix models represents the partition function of 
2d gravity, cf. eq. (1.2.6)). As A ^ Ac, we see that A ^ 00 so that we may rescale 
the area of an individual triangle e to zero, thus giving a continuum surface with 
finite area. Intuitively this is clear, by tuning the coupling constant to the point 
where the perturbation series diverges the sum becomes dominated by diagrams 
with infinite number of vertices, which is precisely what we need for the continuum 
limit. 

The above reasoning constitutes by no way a proof that the matrix model in 
the continuum limit coincides with continuum 2d-gravity. However, one is able to 
compare properties of the partition function and correlation functions calculated 
by matrix models with results directly computable in the continuum formulation of 
Liouville theory. This gives implicit confirmation that the matrix model approach 
is sensible. 

One such quantity computable in Liouville theory is the scaling behaviour of the 
Volume of Sg^^i of eq. (1.2.2) 



Vol(I]g,A) = e^^^A(i-f) Wg, 



(1.2.9) 



where ^c and Wh are undetermined constants. Here the critical exponent 7^4,. is 
called the string susceptibility. The above formula even holds true if one couples 
conformal field theories with central charge c to Euclidean gravity. The Liouville 
theory prediction for the string susceptibility then is [11] 

= 3^(c-l- V(c-l)(c-25)). (1.2.10) 

By plugging eq. (1.2.9) into eq. (1.2.2) and integrating (m^ the area A we find ^. 

,{2-2g)/G _ ^J(E}- ) (1-9) r[ {g - 1)^ 



3=0 



tr-%...)]Wg. (1.2.11) 



Thus we see that the scaling behaviours of the matrix model in eq. (1.2.7) and that 
of quantum gravity in eq. (1.2.11) are identical with 70 — = —1/2. Comparing 
this with eq. (1.2.10), shows that we in fact are describing the c = theory, i.e. 
"pure" 2d-Euclidearrglavity. This result may be further confirmed by comparing 
planar correlation func ions. 

^Simply use the definition of the gamma function r{x) = dt t^~^ . 
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Double Scaling Limit 

There is yet another Umit we have to discuss. As a matter of fact the matrix model 
can only be solved in the iV — > cx) limit, as we shall see in the next chapters. But 
then due to eq. (1.1.16) only the planar contribution to the free energy survives. 
However, the successive contributions Fg all diverge at the critical value for the 
couplings X — Xc- This suggests that if we take the limits iV — > oo and A ^ Ac not 
independently, but together in a correlated manner, we may compensate the large 
N high genus suppression with a A — > Ac enhancement. This results in a coherent 
contribution from all genus surfaces. 

We already saw the leading singular behaviour of Fg in eq. (1.2.7), i.e. I 

Fg[X] = fgiX-Xc)^'-''"^^'-<^\ (1.2.12) 

Then in terms of 

a = N-^X ~ Xc^"-^ (1.2.13) 
the genus expansion (1.1.16) may be rewritten as 

oo 

N'F[X]^Y.''''' fa- 

9=0 

The double scaling [13] limit is thus obtained by taking the limits — > oo, A ^ Ac 
while keeping fixed the "renormalized" string coupling constant a of eq. (1.2.13). 

The above limit may be performed in the case of general potentials as well. In this 
case the additional degrees of freedom can be used to fine tune the couplings in such 
a manner as to adjust an alternative string susceptibility, i.e. at the m'th "multicrit- 
ical" point we find 70 = — 1/m. It may be shown that these theories then correspond 
to the (2, 2m — 1) non-unitary, minimal conformal field theories ^ coupled to quan- 
tum gravity and having the central charge c = 1 — 3 (3 — 2m)^/ (2m — 1). 



3. Virasoro Constraints and the Loop Equation 



matrix model. The 
at the integrabil- 



We now derive a set of constraint equations for the Hermitian 
generators of these constraints obey the Virasoro algebra, alludibg- 
ity of the model as well as a conformal field theory formulation of it. This structure 
sets the basis for a supersymmetric generalization of the Hermitian one matrix 
model to the supereigenvalue model. 



^These are classified by the pair (p,q) with c = 1 — 6(p — q)-^ /pq. 



□ 



□ 



1.3. Virasoro Constraints and the Loop Equation 

Consider the eigenvalue model (1.1.4) under the shift of integration variables 
Xi Xi + e Xi"'^^ , with e infinitesimal and n > —1 free. The resulting expression 
proportional to e must equal zero. One finds 

n 

6 ( iV ^ fc .g, ^ A,^+" + Y^iY^ X--') (^ A/) ) = n> -1. (1.3.1) 



fc>0 



fc=0 i 



Of course this analysis may also be performed in the matrix formulation of eq. 
(1.1.1). 



To see this look at the following variations 
5 



and 



For the van der Monde determinant we have 5 A(A) = A{A) 



and thus 



iy^j ^ fc— i j i 

Putting these three equations together yields eq. (1.3.1). 

The Schwinger-Dyson equation (1.3.1) may be recast in the form [14] 

L,^Zn^O n>-l, (1.3.2) 
where is a differential operator in the coupling constants 



1 " 



k>0 



(1.3.3) 



In the case h 
of a closed 



~1 the second sum simply vanishes. The operators iL„ are generators 
)set of the y ir^soro algebra 



[/C„j , £jn ] — (m — n) Ljn+ri 



(1.3.4) 



without central extension (remember that n,m> —1). This is the reason why 
(1.3.2) are called the Virasoro constraints, in bosonic string theory physical states 
obey an analogous set of constraints. 



I: The Hermitian Matrix Model 



Conformal Field Theory Formulation 

Given a complete set of constraints on a partition function which form a closed 
algebra, one might now ask the inverse question: What is the integral representation 
of the partition function obeying these constraints? In the case of the Virasoro 
algebra it is natural to look for an answer to this question in form of a correlation 
function in a conformal field theory [15] . The methods introduced in this subchapter 
can be easily generalized to the case of the super- Virasoro algebra considered in 
chapter III and lead to the construction of the supereigenvalue model. 

Just as the generators of eq. (1.3.2) the modes T„ of the energy-momentum 
tensor of a free scalar field obey the Virasoro algebra. Let us thus consider the 
simplest possible conformal field theory, a holomorphic scalar field 

cj,{z)^q+p\nz + Y,'-^^\ (1-3.5) 
with the commutation relations 

[Jn,Jm] = nSn+m,0, [ 9 , P ] = 1- (1.3.6) 

Define the vacuum states 

Jfe |0) =0 (A^l J_fc = fc > 

p|0) = 0, (1.3.7) 

where |A^) = exp[v^iV g] |0) and q is anti-Hermitian. The energy-momentum ten- 
sor is given by 

T(z) = i :[a0(z)p:=^r„z— ^ 

nez 



□ 



Tri = ^ J-k Jk+n + ^ JaJb U > 0, (1.3.8 



2 

A:>0 a-|-6— n 



and we define a Haniiltonian by 



H{g,) ^J-J2g,Jk = ^<f ^ Viz) dcbiz), (1.3.9) 



where V{z) = J2k>o9i^ '^^ ■ With these definitions one shows that 



(iV| exp[i?(5fc)] 



- {N\ exp[i?(5fc)]T„ 



(1.3.10) 



1.3. Virasoro Coi 



straints and the Loop Equation 



So any operator S satisfying 



[T„,g]=0, n>-l, 



will lead to ^ 



L„ {N\ exp[H{gk)]9 



m 



= 0, 



(1.3.11) 



(1.3.12) 



as Tk |0) = 0. Therefore any nonvanishing correlator of this form is a candidate for 
the partition function 

ZN = {N\cxp[Higk)]9\0), ril.3.13) 

obeying the Virasoro constraints. Finding the operators S satisfying eq. (1.3.11) is 
an internal problem of conformal field theory. The solution is given by an arbitrary 
function of the screening charges Q± which are of conformal dimension zero 



Q±= j) duj : exp[±\/2(j){u;)] 



(1.3.14) 



Choose S = Q+ in order to get a nonvanishing correlator in eq. (1.3.13) ^. We then 
have 

Zn — {N\ : exp [-= i V{z)d(j){z)]: Yl f dz, : exp[V2 (j){zi)] : |0). (1.3.15) 
v2 J Co 

Using the operator product expansion 4>{z) 4>{z) ^ ln(z — z) this is evaluated to 

N 



(1.3.16) 



After deforming the contour integrals to integrals on the real line we thus recover 
the Hermitian matrix model in the eigenvalue representation of eq. (1.1.4). 

To obtain cq. (1.3.16) from eq. (1.3.15) use : : : := : e^+B : , therefore 

JV „ „ AT 

^N=(j[j) dz,'^ {N\ -.eM^ J) V{z)d^{z)+V2,t>{zi)]:Y[:cMV2,t>{z,)]:\0) 



iV „ . N 

Y^j) dz,{N\:eM^f V{z)d<P{z) +V2{^{zl) + ,|>{z2)}]■.'[[■.cM^H^^)]■■ \0) 
i^l-JCi J Co 



*Note that we have extended the definition of T„ to n = — 1 by dropping the second term in 
eq. (1.3.i 



^More generally one could take Q^"*"*^ S*^, for more details see [15]. 



I: The Hermitian Matrix Model 



iV „ „ 

= T\f dz, {N\ -.expl^ (b +v/2 V0(^,)] : |0> 

i i<j 

N 

= Y[f dz,{N\ :oxp[v^Arg]: |0> exp[ ^ ] - 

i = l '^i i i<j 

and as (A''| : C'Kp[\/2N q ] : |0) = 1 we have proven eq. (1.3.16). 



Loop Equations 

An immediate consequence of the Virasoro constraints is the loop equation, an 
integral equation for the one-loop correlator W{p) of eq. (1.1.6). Multiply eq. (1.3.1) 
for n with 1 /p""*"^ and sum up the obtained equations. The resulting geometric series 
^ may be performed to give 

i i 

Note that generally (E, OiX.f )conn. = ( E^ 0{X^f ) - ( OiX,) f and hence 

{^^y^]^)^W{p,p) + N^W{p)\ (1.3.18) 

with the definition of the loop-correlators of eq. (1.1.8). In order to rewrite the first 
term of eq. (1.3.17) in an integral form introduce the eigenvalue density 

'^W = ^E(^(^-^0>- (1-3.19) 

i 

The eigenvalue density has the remarkable property that it vanishes outside a sup- 
port of [y, x] on the real line in the limit of infinite matrix size, N — > cxi. This can 
be seen by studying the saddlepoint evaluation of the eigenvalue model (1.1.4) in 

''The subtlety of p larger or smaller than A; may be removed by directly considering a shift 



1.4. The Solution 



this limit [10]. It is precisely this limit which we consider in the foiipwing, keeping 
all subleading terms in of the genus expansion. In this sense we will be able 

to solve the matrix model pcrturlja|jvely by expanding around tbeJpoint N = oo. 
Each term of this expansion will be some function of rtinc couplings gk, which can 
be determined exactly. So the Hermitian matrix modelwill turn out to be solvable 
nonperturbatively in couplings but perturbatively in 1/N'^. 



By using the eigenvalue d 
comes 



msity of eq. (1.3.19) the first term of eq. (1.3.17) be- 



^ P-K' J p 

= N^ ( d\p{\) 



-A 

duj 1 V'{uj) 
2'Ki uj — X p — uj 



(1.3.20) 



»ver p lies 



Here the curve C of the contour integral runs around the cut [a;, y]. JV for 
outside the curve. The integral over A may now be performed in eq. (1.3.20) yielding 
the one-loop correlator W{Ld). Combining this with the result of eq. (1.3.18) finally 
leads us to the loop equation of the Hermitian one-matrix model [16, 17, 18] 



2lT p — LU 



W{pf = ^W{p,p) 



(1.3.21) 



This integral equation basically captures the whole Hermitian matrix model, 
eq. (1.1.11) we have W{p,p) 



Due 



5/5V{p)W{p). So the loop equation (1.3.21) is 
closed equation for the one-loop correlator TmXpI, determining this quantity. As 
explained in section 1 knowing W{p) then means knowing all correlators of the 
model. 



4. The Solution 



There are various methods to solve the Hermitian one matrix model. Two of the 
most prominent ones are the saddlepoint evaluation [10] and the method of orthog- 
onal polynomials [19]. The most effective solution, however, is based on an iterative 
procedure to solve the loop equation (1.3.21) genus by genus [8]. This method is 
the one which has found a generalization in the supereigenvalue model, so let us 
briefly review its basic concepts. i — 

The loop equation (1.3.21) may be solved by mfkiiig use of the genus expansion 
of W{p) I— 

OO - 

9=0 
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I: The Hermitian Matrix Model 



Plugging this expansion into the loop equation (1.3.21) and comparing terms of 
common order in 1 /N'^ leads to a coupled hierarchy of equations for the Wg {p) . To 
leading order in 1 /N'^ one simply has 



dco V'iio) 

Wo w 

Q Ztt p — uj 



Woipf 



□ 



(1.4.2) 



From the definition of W{p) in eq. (1.1.6) we know that asymptotically 
Wo{p) = l/p + 0{p~^) for p — + oo. If one additionally assumes that the singularities 
of Wo{p) consist of only one cut on the real axis ^ one finds [16] 



^^^^^ j; dw V'ico) iip~x){p~y)i 1/2 



2tt p ~- oj 



(1.4.3) 



where the endpoints of the eigenvalue distribution x and y are determined through 
the matrix potential in the following way: 



duj 



C 27r ^{u; -x){lj- y) 

duj u)V'{uj) 
C 27r y/{uj -x){uj- y) 



0, 



= 2, 



(1.4.4) 

□ 

(1.4.5) 



which are a direct consequence oi W [p) = 1 / p + 0{i ^) 



Let us verify the planar solution (1.4.3): 



dw 



(P - a;) (p - y) 



dz_ V'{u)) V'(z) 

2n {p~uj)(p~z) [(i^ - x){u) - y){z - x){z - y)Y/'^ 



Now rewrite {j> — x) = {p — u)) + (w — x) in the numerator. By eq. (1.4.4) only the (a; — x) terms 
survives. Doing the same for (p — y) gives 



WoipY 

Using 



2^ 



dz V'{uj)V'{z) {ui-x){z-y) 
2tt {p-uj){p~z) [{u; - x){ui - y){z ~ x){z ~ y)]^/^' 



(1.4.6) 



1 



1 



1 



1 



1 



p — UJ p — Z UJ — z p — UJ P — z 
in eq. (1.4.6) and renaming a; *-> z in the resulting second term one finds 



Wo{pf 



duj 
2^ 



dz V'{uj)V'{z) {uj-x){z-y) + {uj-y){z-x) 
2tt {p -uj){uj- z) [{uj- x) (uj - y) {z - x) {z - s/) ' 



^This follows from the relation of Wo(p) to the eigenvalue density of eq. 
p(z) = lim.^o 2^ [Wq{z - ie) -Wo(z + it)]. 



(1.3.19): 
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1.4. The Solution 



(The process 'of pulling the curve C2 over Ci gives no contribution). Now we do the same trick 
again and write (z — y) = (z — uj) + {w — y) as well as (z — x) = (z — uj) + (w — x) to get 



Woipf 



duj 
2^ 



C2 



dz_ V'{uj) V'{z) 

27T [p — Uj) {u! — z) 



(uJ -x){uj- y) 



{z -x){z- y) 



1/2 



duj V'(uj) 
2-K p — UJ 



W{uj), 



which proves eq. (1.4.2). 



An interesting consequence of the plamr solution (1.4.3) is the rather universal 
form of the planar two-loop correlator WTj{p, p) 



Woip.p) 



[x - yf 



16 (p — a;)2 {jp — y)2 



(1.4.7) 



depending on the matrix potential only through the endpoints of the (iig<:nvalue 
distribution. 

Let us now turn to the higher genus contributions of W{p). By plugging the 
genus expansion (1.4.1) into the loop equation (1.3.21) it appears that the Wg{p) 
for 5 > 1 obey the equation [8] 



3-1 



V' O Wg{p)=Y,Wg^{p)W^{p) 



□ 



5V{p) 



Wg-l[p), 



(1.4.8) 



where we have introduced the linear operator V by 



V o /(p) 



duj V'{uj) 
2t: p ^ UJ 



f{uj)^2Wo{p)f{p). 



(1.4.9) 



In eq. (1.4.8) Wg{p) is expressed entirely in terms of the Wgi , g' < g. This makes it 
possible to develop an iterative procedure to solve for Wg{p) as done by Ambj0rn, 
Chekhov, Kristjansen and Makemko [8]. One basically has to find a way to invert 
the operator V of the left hand sldaof eq. (1.4.8). It turns out that this is effectively 
done by expanding the Wg (p) in a set of basis functions of V augmented by a change 
of variables from coupling constants to moments of the potential. As we shall 
encounter precisely the same problem in our discussion of the iterative solution of 
the superloop equations for the supereigenvalue model in chapter III, let us postpone 
the analysis of this point. 



□ 



Chapter II: 

The c = -2 Matrix Model 



To construct supersymmetric versions of matrix models it is necessary to augment 
the one matrix model of chapter I by fermionic degrees of freedom. There are various 
ways to do this. On the level of matrices one could think of replacing the Hermitian 
matrices M by supermatrices M of type {N\M) 



i,j = l,...,N 



^^^=yC yZ) a,P^N + C--\N + M, (2.1) 



where the Xij and Yap are Grassmann even and '^ia and '^ai are Grassmann odd 
quantities. A supermatrix model is then naturally defined by the partition function 
of eq. (1.1.1) where one replaces the trace by a supertrace. This model has been 
studied in refs. [20], but the outcome is rather disappointing: The partition function 
Zs{N\M) of such a supermatrix model of type {N\M) turns out to be proportional 
to the partition function Zg{N — M|0), which is nothing but an ordinary matrix 
model built out of {N - M) x {N - M) matrices ^ 

In this chapter we shall address an alternative approach to this problem. Guided 
by the superspace formulation of supersymmetric field theories, one augments the 
zero dimensional space of the Hermitian matrix model by two anticommuting co- 
ordinates. The fields living in this space will be matrix valued superfields, whose 
component expansion consists of two bosonic and two fermionic matrices. Think- 
ing of the interpretation of anticommuting coordinates as negative dimensions, this 
model should be considered as a d = —2 matrix model. The planar triangulations 
of closed surfaces in c? = — 2 dimensions was first studied on combinatorial grounds 
by Kazakov, Kostov and Mehta [21]. The matrix model approach to this problem 
was introduced by David [22] and further studied in refs. [23, 24]. And in fact in 
its scaling limit the matrix model does describe the coupling of c = —2 matter to 
2d gravity. 



A negative {N — M) corresponds to a change of sign in the exponent of eq. (1.1.1). 



In our analysis of the model we shall stress its supersymmetric structure and 
discuss the emergence of a Nicolai-map to a Gaussian d = Hermitian matrix 
model. The diagrammatic interpretation of the model is investigated, leading us to 
its equivalence to random surfaces decorated by fermion loops. The one-, two- and 
three-point functions are computed in the case of a quartic potential, which has not 
appeared in literature so far. In the critical regime of this solution we recover the 
scaling exponents of c = — 2 Liouville theory. Finally as a nice application of the 



solution we present the counting of the first few supersymmetric graphs, si nilar in 
fashion to the counting of bosonic planar diagrams in ref. [10]. 



1. The Model 



□ 
□ 



In chapter I we studied Hermitian matrices living in zero dimensional space. As a 
method of supersymmetrization of matrix models, let us consider matrices living 
in superspace. In order to maintain exact solvability we introduce two fermionic 
coordinates 9 and 6, but keep thte-^ro dimensional bosonic space. A matrix valued 
superfield O will then have the following component expansion 



<S< = ^ + ^( 



(2.1.1) 



where # and F are Hermitian N x N matrices and the >P' and <P' are N x N matrices 
with complex Grassmann odd entries. We have = 't' ajndj 9* = 9. Note that O 
is Hermitian. 



The matrix model built out of O is given by the partition 



Fun ction 



Zcj>[gk] = I 'DO cxp 



-NTiS[(t>] 



(2.1.2) 



As measure we take V^VFV^'D^ , where the measure of the Hermitian 

matrices $ and F are as m eq. (1.1.2) and for the fermionic matrices we take 



N 



n 



af3 d'l'pa 



(2.1.3) 



a,/3=l 



II: The c = -2 Matrix Model 

The action of the O-matrix model of eq. (2.1.2) reads 



2.1. The Model 



Supersymmetry Transformations 



k=0 



After integrating out the auxihary matrix F we find that the effective action of eq. 
(2-1-4) (2.1.7) is invariant under the following set of (global) supersymmetry transforma- 
tions — 



with the "superspace" derivatives D = d/dO and T) = d/dO. Bearing in mind the 
interpretation of anticommuting coordinates as negative dimensions it should be 
clear why eq. (2.1.2) is considered as a d = — 2 matrix model. 

The obvious extension of this model to one bosonic dimension, i.e. considering 
component fields ^(x), ^{x)^ ^{x) and F(x) depending on oie variable x, is known 
as the Marinari-Parisi superstring [25] ^ . 

After performing the integrals over 6 and 9 in the action (2.1.4) and using cyclicity 
under the trace one finds 



1 



a/3 



■eV^'(*)a/3 



1 



where e and e are anti~commuting parameters. 



(2.1.9) 



TrS'[cD] = -TrF^ + Try'(#)F + ^/cgfe ^ Tr #° •P' 

Q + fc=fe-2 



This is easily verified. Consider the variation of eq. (2.1.7) under the transformations (2.1.9), 
(2.1.5) then 



k=0 



using the general matrix potential 



cancels with 



3=0 



-^Tte'l'V'i'P) V"(#) - ^TrV'i^) V" (^) 'P e + {S^'s} . 



a + 6=n-2 



As the auxiliary matrix F enters in the action Tr S'[ O ] of eq. (2.1.5) only quadrati- 
cally, a shift in integration variables Fap i-^ F^p + 1/2 V'{^)ai3 lets the integral over 
F decouple. The Jacobian associated with this shift ipTIhity, hence the integral over 
the auxiliary matrix F can be performed yielding an |y_iependent constant cat. We 
are thus led to the effective action 



1 



(2 



and the partition function now has the form 

JNxN 



riceiT-^I 



(2.1. 



The next obvious thing to do is to integrate out the fermionic matriceST-^Let us, 
however, first investigate the symmetries of this effective action. 



except for the variations {(5^'s}, which cancel by themselves. See this by looking at the part of 
the transformations (2.1.9) proportional to e 

{5#'s}e= ^ Tr<5f(#'")'J'<f''lf' + Tr#''^<5f(#'')«J>- 

a-\-b-\-c—n — 3 

the calculation for the part proportional to e works analogously. 

If one goes back to the form of the action including the auxiliary matrix F-Qf eq. 
(2.1.5), the supersymmetry transformations become linear in the fields: 

S^af} = -eFap 
= -e Fap 

SFc,f3 = 0. (2.1.10) 



^The question whether the continuum limit of this model describes a theory with target space 
supersymmetry is still unsettled [26]. 
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II: The c = 



-2 Matrix Model 



□ 2.1. The Model 



Note that the measure PlP' is invariant under the supersymmetry trans- 

formation (2.1.9). This immediately follows from the form of the Jacobian of the 
transformations (2.1.9) which is zero in linear order of e and e. 

Integrating out the Fermions 

Just as in the Hermitian matrix model the action TrS'^g of eq. (2.1.7) is in- 
variant under the simultaneous U{N) transformations # i-^ U^^U, ^ i-^ W^U 
and «P' C/t«P'C/. This may be employed to diagonalize the Hermitian matrix 
# = [/t ^j-, If, with #o = diag {(pi, . . . , (Pn)- We again pick up the van der Monde 
determinant nQ</3('/'a — fp)'^ the Jacobian of this transformation. As the trans- 
formation matrix U does not depend on the fermionic matrices ^ and <P', the Jaco- 
bian of their transformation is unity. We thus obtain the semi-eigenvalue picture 
of the O-model 



[ g/c ] 



N 



a=l a, (3=1 a<f3 



NS,s[gk]\, (2.1.11) 



where the constant cat now contains the integral over U{N) and the constant cat of 
eq. (2.1.8). The partially diagonalized effective action reads 

Se{i[9k] = lj2^'i'Pc.f +Y.kg, J2^a^'-P^f^^P'^- (2.1.12) 

a k a+b=k—2 a(3 

In this form the integration over the fermions may be performed rather easily. 
To do this introduce the abbreviation 



vi4 



a+b=k-2 



if a = /3 




(2.1.13) 



We then have 

N 



/ Yl d^c^P d^c^P exp [- Cap ^ccP ^I3a 

a, 0=1 a,l3 



N 



n[C"« 



a.p 



a<l3 



<fa - fp 



(2.1.14) 



This is a rather astonishing result, as the van der Monde determinant will now 
drop out of the partition function. Plugging eq. (2.1.14) into eq. (2.1.11) yields the 
purely bosonic partition function 



-o[gk\ 



- CN 



1 2 



exp 



N 



So by making the N independent substitutions 

K^V'{ipa) => dXa ^V"{ipa)d(pa a = 1, 

the partition function [ ] takes the Gaussian form 



-ttifffcj = Cat 



N 

n 



V'(oo) 

dxA \\{\o 



EN 



(2.1.15) 



(2.1.16) 



(2.1.17) 



note that this is true for any continuous function V'{Lpa). 

We conclude that in the case of unbounded potentials V'{—oo) — V'{oo) the 
partition function Zc^ vanishes. On the other hand for potentials with a lower 
bound, i.e. V {—oo) = —oo and V'{oo) — oo, the O-matrix model is proportional 
to a pure Gaussian Hermitian matrix model 



2o[gk_ 



CN 



VA exp 



N 



■TiA-" 



r,N ^N' 



(2.1.18) 



and is completely independent of the coupling constants g^! The result 2^ tt^^ is 
obtained by collecting all N dependent constants and making use of the result for 
the Gaussian partition function [19]. 

But this means that the partition function can obviously not have any critical 
behaviour. Moreover any integrated supersymmctric expectation value vanishes 



Jl d^c^p d'^ap Y[[^+ Cap '^ap ^pa 

a,p=l a,P 



d0 d^TrO" 



0, 



(2.1.19) 
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II: The c = 



-2 Matrix Model 



2.2. Feynman Diagrams 



which is a direct consequence of differentiating Zt^ of cq. (2.1.2) with respect to 
(7„. On the level of the effective action (2.1.7) this Ward identity reads 



a-\-b— 71 — 2 



(2.1.20) 



Nevertheless non-supersymmetric correlation functions like ( Tr 0*^ ) are non-trivial 
and display critical behaviour as we shall see in the subsequent sections. 

Nicolai Map 



I I 



>x 



N g 



2 


: > 


N 




N 




— 8 
2 


VaV 



1 

N 

N 2 
— 8 
4 



A', 



X 



L 



Ng 



Index Loop: 



Fermion Loop: 



The observed situation is very reminiscent to the Nicolai-map [27] of globally su- 
persymmetric field theories: Integrating out the fermions in a supersymmetric field 
theory yields a determinant and an effective bosonic action. The vanishing of the 
vacuum energy of these theories alludes at their equivalence to a free theory. As a 
matter of fact one can show that there always exists a map from a free, Gaussian 
action to the effective bosonic action of the supersymmetric theory, whose Jacobian 
is identical to the determinant of the fermionic integration. On the level of correla- 

tion functions this implies that an operator expectaticpryalue of the supersymmetridZH To be specific we consider a quartic potential of the form 

cated operator expectation 



Figure 3: The Feynman rules for the quartic O-matrix model. 



theory may be calculated as a generically rather compl: 
value in a simple free theory. 

This scenario directly translates to the O-matrix model. The correlator 
( Tr ) may be evaluated in the Hermitian matrix model with Gaussian measure 
Sq — — A^/4Tr by using the inverse map of eq. (2.1.16) 



1/((D) 



.9^ 



(2.2.1) 



Tr**^ ) = 



Tt[V'-\A)]''\ 



where 



denotes the expectation value in the free 



(2.1.21) 



troory 



OiA))) = ^[ 2?AO(A) exp[-^TrA2]. (2.1.22) 

We shall exploit this relation in sections 3 and 4 for the calculation of bosonic one-, 
two- and three-point correlators. 



Note that the quartic poteni ial is the simplest non-trivial bounded potential to 
consider, i.e. satisfying y (— ob)J= — oo and V'{oo) = oo for g positive. To our mind 
this point has not found adequate attention in the literature [22, 23, 24] where a 
cubic potential is assumed and the boundaries of intefrStion in eq. (2.1.17) are set 
by hand to — oo and oo. 

With this quartic potential we are led to the effectiEeHaction for the matrices ^, 

and via eq. (2.1.7) 



Tr5eff, = iTr*^ 



igTr#4 
2'^ 



- g2 Tr 



+Tr<l^<i^ + 5|Tr#2<P'!P' + Tr*«P'#!P' + Tr«P'#2^|. (2.2.2) 



2. Feynman Diagrams 

Let us derive the Feynman rules for the O-matrix model. This will lead to a 
geometrical interpretation in terms of "fat" graphs, as we saw in the Hermitian 
matrix model in chapter I. 



The bosonic and fermionic propagators read 



[ ^a(3 ^-yS )g=0 = ^aS 



J7J 



[ ^aP ^^75 )g=0 = ^Q<5 SjP^, 



(2.2.3) 



and again they are represented by double lines each one corresponding to the sep- 
arate propagation of matrix indices. Moreover one has to assign a direction of 



15 



II: The c ^ 



-2 Matrix Model 



2.2. Feynman Diagrams 



propagation to the fermions due to a necessary canonical ordering of the matrices 
^ and ^. 



The form of the propagators in eq. (2.2.3) is directly derived by 



2 



^7 



oxp 



N 



N 



4 Op J 



and 



5^/ 



exp 



TV 



via partial integration. 

With the interactions of eq. (2.2.2) there will be a bosonic four-point vertex 
contributing {—N/2 g), a bosonic six-point vertex contributing {—N/ 4 g^) and four- 
point Yukawa vertices contributing {~N g). Again each loop of internal index will 
yield a factor of N — 6aa- Moreover there is a factor of (—1) for every closed 
fermionic propagator loop. To see this consider a fermionic loop in an arbitrary 
diagram built out of n Yukawa vertices. We then have to evaluate the following 
expectation value, where the dots stand for bosonic propagators running out of the 
Yukawa vertices and all matrix indices are suppressed 




Figure 4: A dual graph of the O-matrix model. 

function as coming from a Hermitian matrix model with a (4, 6) vertex inter- 
action, where one draws / loops running over np 4-point vertices on the diagrams 
and weights each resulting diagram by a factor of (—1)'. (This is true as np boson 
propagators are turned into fermion propagators yielding a factor of 2"*', which 
is exactly the factor needed to convert np bosonic 4-point vertices into Yukawa- 
vertices). Note that the loops drawn on the graph are not allowed to intersect or 
touch. Hence 



/ Ok 
5 ^ 



^2 -^2 



5^4 



5fn 



(-){.. 



^2 -^2 



One thus gets a factor of (—1) for every closed fermionic propagator loop. 

Hence a diagram consisting of rt4 4-point and rig 6-point bosonic vertices, np 
Yukawa- vertices, Eb boson propagators, Ep fermion propagators, L index loops 
and I fermion propagator loops contributes the factor 



Diag 



(_)' (-|)"' (-y)"' (-5)"'' iV"^+"«+"'=- (^^Y""^^^ 2^« (2.2.4) 



(0 ( 

herm^^'^ 



l/4,54-5/2,56-5V4) 



(2.2.5) 



1=0 



in the notation of cq. (1.1.3). Interestingly enough the loop drawing lets the sum 
in eq. (2.2.5) trivialize. On the level of the dual graphs we see that each diagram 
of represents a discretized surface built from squares and hexagons with loops 
painted across squares. 



to the partition function. Using the Euler relation one again has the overall genus 
dependent N'^^'^s weight for a graph of genus g. One may envisage the partition 



As a nice check of the triviality of the partition function Z,f, of eq. (2.1.18) let us evaluate the 
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0{g) vacuum graphs: 




[1] 



-^9 



2( 



[2] -^3 



[1] -^9 



n) (n) 

(tt) (w 



iV3 = 4 c, 



N = 2g 



Keep in mind that the six-point vertex is of order g^. The numbers in [. . .] denote the combinatorial 
factors. We see that the sum of all vacuum diagrams to first order g vanishes, as they should. 
Note the genus controlling factor of A'^^ in the above graphs. 



3. Bosonic One— Point Functions 

As mentioned in section 1 bosonic one-point functions of the form ( Tr ) might 
display critical behaviour despite of the triviality of the partition function . 
To calculate these quantities one makes use of the Nicolai-map of eq. (2.1.21) to 
map these correlators onto correlators of a Hermitian matrix model with Gaussian 
potential. 

Let us again assume the quartic matrix potential of cqs. (2.2.1) and (2.2.2). As 
it is symmetric all odd bosonic one-point correlators (Tr#^"+^ ) vanish. For the 
even correlators one has to solve the equation 



(2.3.1) 



for (p{X) — V'^^(X). For g > there is only one real solution to eq. (2.3.1), whereas 
for 5 < we have three. Picking the branch selected by smooth continuation of 
g —g we find 

v'-\x) = 



x + 



A2 



s 

27 



^A2 + A 
4 27 



(2.3.2) 



where we have introduced s — sign(5). 



2.3. Bosonic One-Point Functions 

As all the correlators (( Tr A^'' )) in the free theory are computable, all one has to 
do in order to calculate the expectation values (Tr^^" ) is to take the 27i'th power 
of the result (2.3.2) and expand the resulting expression in A. The outcome of this 
straightforward calculation is 

n / r, \ ^ ( ^ 1 



k=l 



(-35J 



(2i + 2)! 

which is true for positive and negative g. 

In order to obtain this result, introduce the abbreviations 



Tryl 



2i+2 ' 



(2.3.3) 



(f) 



21 V 4 , — ,27 



Then y'-i(A)2" = [(+) + (-)V"/\g\''' 

n 2n 



Ma2 + -. 

4 27 



(+) + (-) 



J2{n + k) ( 3) 



:femg (+) ■ {—) = —s/3 we have 

n — k 



+ (-) 



2fc 



2n 



(2.3.4) 



Consider the case s = 1. With the expansion [28] 



{x + ^/l + x^)^^ + (x- ^/l + x^)1 = 2(1 + ^ 

i=0 

and eq. (2.3.4) one arrives at 



g^[g^-22]...[g2-(2»)^] ^2,+2^ 
(2i + 2)! ' 



n 

y'-(A)-=;^2(-)'-(„2;^) 

With the help of the identities 



(39)-"+^m[fc2_(302] 



(39)' 



(2i + 2)! 



(2.3.5) 



and 



!:(..")'-''-;(») 

n 
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which one may prove, one recovers the result of eq. (2.3.3) from eq. (2.3.5). A similar analysis is 
possible for s = —1 yielding the same result as eq. (2.3.3). 



Biatrix model 
ef. f8l sketched 



Gaussian One Point Functions 

In order to calculate the one-point functions of the Gaussian 
(( Tr A^'' )) one can use the iterative solution of the loop equations of ile 
in chapter I. However, due to the simplicity of the Gaussian model itis more con- 
venient to employ the method of orthogonal polynomials [19]. Here the evaluation 
of ((Tryl^'^)) for all genera reduces to an integral involving Hermite polynomials. 
This problem was solved by Kostov and Mehta [23] and we simply cite their result. 
For the matrix potential 1/4 yl^ of eq. (2.1.22) the result reads 



1 

N 



2i+2 



yi+1 



(i + 1)1 {1 + 2)1 



where J'„(2iV) is a polynomial in (2iV)" 



[„/2] 



y„(27V) = ^ 



(27V) 2J ' 



whose coefficients Unj are defined by the recursion rdlallion 



JiaJic 



Qn+l, j 



k=2j-l 



k{k + 1) ak-i,j-i, 



(2.3.6) 

(2.3.7) 

□ 

(2.3.8) 



with a„o = 1- Note the genus expansion recovered in ec~\l.3.7). The fact that one 
has a recursive solution for higher genera contributions-ae not very astounding if 



one recalls the general iterative solution sketched in chapter I. 



a„o = 1 
1 
3 

1 

90 
1 

5670 
1 



a,nl '■ 



an2 ■ 



a„4 : 



(n + l)n (n ~ 1) 
(n + l)n{n- 1) (n - 2) (n - 3) [5n - 2] 
(ra + 1) ra . . . (n - 5) [35n2 - 77n + 12] 
(n + l)n 



340200 



(n - 7) [175n^ - 945n2 + 1094n - 72] 



: (n + 1) n ... (n — 2j + 1) [ Polynomial of degree ] 



(2.3.9) 



Results for all Genera 

Putting the results of the previous subsections together, i.e. plugging eq. (2.3.6) 
into eq. (2.3.3), yields the bosonic one-point correlators of the O-matrix model 
with quartic potential for all genera 



G-n 



N 



Tr^ 



2n 



oo ^ 



h=0 



]^2h ■ 



Let us define 



k=l 



2n 

n — k 



j+n-l 

n [(30= 

1=0 



(2.3.10) 



(2.3.11) 



Note that Aj simplifies to = 3 ^ (3j + l)!/(j!). Then we have the genus h cor- 
relators G'Jj(g): 



2i+n+l 



{j + n)l {j + n + 



j^K(7l,h) 
OO 

G'i= E ^1 

j=K,(n,h) 
oo 

Gl= E ^1 

oo 

Gt= E ^" 

j=K{nJi) 



2 j+n-l 



(-3.9)^ 



1 



o (-35)^ 



^ U + n)lij + n-2y. 3 

2J+"-3 [5{j + n)-2] 



(-3.9)^ 



(j + n)l {j + ji - Ay. 90 

2J+"-5 [35(j + n)2-77(j + n) + 12] 



U + n)\{j + n-6)\ 

2 j+n-7 

{j +n)\{j + n~8y. 



5670 



(-35)^ 



;i75(j -f n)^ - 945(j -t- n)^ -|- 1094(j + n) - 72] 



340200 



(-35)^ 



Gh 



oo 



A'- 



2j+n-(2h-l) 



. Vm''' u + nyij + n + iy 



aj+nJi h^9)\ 



(2.3.12) 
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where we have introduced k = max (0, 2h — n) as the lower bound of summations. 
The correlators all have the same radius of convergence of jgd — 1/54. Moreover 
the planar G'^(g — gc) are seen to converge at the critical point, whereas the higher 
genus correlators G^^{g = gc) diverge at the critical value of the coupling constant. 

Let us now compute the asymptotic form of the one-point correlators G'^{g). The 
factors A" of eq. (2.3.11) behave as 



(j+n-iy. \ (2.3.13) 



for J 3> 1. Using this we compute the scaling behaviour of the one-point functions 

2j+n-2h+l 

U + ny. {j + n - 2/i)! 



_2«-2'«+i9«-i (j + n ^ 1) (j + n - 2) . . . [j + n-2h + l) f^_^ ^ , 

with gc = 1/54. Hence the correlators scale like 

G:i(ff)~(5c"<?)'-^^ (2.3.15) 

and are independent of n. The bosonic one-point correlators Gjj are to be inter- 
preted as sum over random surfaces of genus h with an n-gon hole as a boundary. In 
fact the form of the boundary becomes irrelevant at the critical point of eq. (2.3.15). 
One now argues that the sum over surfaces without boundaries is obtained by in- 
tegrating eq. (2.3.15) over g. The resulting scaling behaviour 

dgG'^{g)^{gc-g)^'~^-'^^^'-''^ (2.3.16) 

may then with the help of eq. (1.2.7) be associated with a string susceptibility of 
7str. = —1, which via eq. (1.2.10) corresponds to a central charge of c = —2. In this 
sense we have verified the claim that the O-matrix model of eq. (2.1 .2) describes 
random surfaces embedded in —2 dimensions. The double seal ng limit of this 
theory and the connection to the Liouville approach of this problem is studied in 
refs. [24]. Let us remark that despite the supersymmetric structure of this model it 
does not show any correspondence to two dimensional supergravity. 

The results for Gjj encode as well all the fermionic one-point correlators 
12a+b=n Tt^'^sP^^'^' if one makes use of the Ward identities presented in eq. 
(2.1.20). 



2.4. Planar Two- and Three-Point Functions 



4. Planar Two— and Three— Point Functions 



The calculation of higher point functions goes along the same lines. Howe ver, in 
this section we shall be less ambitious and study only planar contributions. We 
keep the quartic potential of eq. (2.2.1). I — I 

For the two-point functions there are now two nonvanishing types of correla- 
tors possible: The doubly even powers in (Tr^^'^^^') and the doubly odd ones in 
(Tr#^'^~-'^ ^2i-i^ Pqj, ^j^g computation of the last type we will have to know the 
power expansion of 1/'~^(A)^"+^ in order to apply the Nicolai-map. This is just 

algebra, take the (2n + l)'st power of eq. (2.3.2) and 
tcome of this calculation is 



another exercise in elementc ry 
expand the result in A. TheLoi 



TrV'-'^{A) 



2n+l 



E(-)'(2fc + i) 




[3i-k~2][3i + k-l] 



( 3g)' /i2/+i 
(2^ + 1)! • 



(2.4.1) 



The next ingredient needed are the two point functions of the Gaussian matrix 
model. Here there is no nice closed form known in the literature. The planar result, 
however, may be obtained rather easily from the planar one-loop correlator of eq. 
(1.4.3) through application of the loop insertion operator as discussed in eq. (1.1.11) 
of chapter I. Specializing to the potential 1/4 TrA-^ one finds the connected, planar 
two point functions 



TrA^^TrA-^' )) =2 

conn 



i + l \ i 



Tr A^^-^Tr A^'-i 



)i+l-3 



through an expansion of Wo{p,q) in the loop parameters p and q. From this one 
may now compose the planar bosonic-two point functions 
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(2.4.2) 



Tr^-^^Tr??^ 



Tr*-' 



Tr#^" 



= (( Tir(^-i(A)2'"TrT/'"i(A)2" 

" conn 
\2m— 1 ta/— 1/ A\2n—1 



-1(A) 



(2.4.3) 



by plugging eqs. (2.3.3), (2.4.1) and (2.4.2) into the right hand sides of eq. (2.4.3). 
We do not write out the results explicitly, as they are lengthy and uoil too instructive. 
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For the evaluation of the planar three-point functions take Wo{p, q, r) and expand 
in p, q and r to read off the planar, connected three-point functions of the Gaussian 
model [29] 



TiA^^TrA^J'TrA^'^ 



conn. V * 



2i 



2i 



2j\ f 2k 



(2.4.4) 



And we may again compose the three-point functions of the O-matrix model by 
plugging these results into the appropriate combinations of eqs. (2.3.3) and (2.4.1). 

Needless to mention that these results may be employed for the calculation of 
higher-point amplitudes involving fermions by making use of higher-point Waird" 
identities. These are obtained from the one-point Ward identity of eq. (2.1.20) jy 
appropriate differentiation with respect to gm, e.g. 



- ( Tr ^ 
2 



n-t-m — 2 



i Tr V'i^) ^"'-^ Tr V'{^) ^"-^ 



\ H (TrF'(#)#"-iTr«f"«i>«f''«P' 



a+b=ri-2 



\ ^ / Tr y (#) - 1 Tr 



E 

o + b=™-2 
c+d = Ti-2 



Tr \['^''\p Tr \p tp 



(2.4.5) 



□ 



5. Supersymmetric Graph Counting 

Quite similar to the classical paper on planar diagram counting of Brezin, Itzykson, 
Parisi and Zuber [10] the results of eq. (2.3.12) may be viewed as a solution to the 
combinatorial problem of supersymmetric graph counting with the Feynman rules 
derived in section 2. 

Let us first consider the planar diagrams of this theory. For the counting of all 
diagrams with 2-bosonic legs we have to write down Gi{g) of eq. (2.3.12) in the 
first few orders of g 



G?(ff) 



lim 4 (Tr*^> 2 - 16.g + 280g2 



(2.5.1) 



Th|(3sn 



terms correspond to the planar diagrams 

2 



O 



-32g 

128^2 

-24.g2 

-325^ 

-645^ 
12852 



16g S : -I65 

512g2 C) O 512g2 

-48.g2 -48.92 

-16g2 y -256g2 S : 280g2 

-25652 (;) C) -256g2 



(2.5.2) 



This constitutes an independent combinatorial check of our results. 

Similarly one can look at the amplitudes with 4-bosonic legs of planar topology 
by expanding G'2(g) 

G°(5)= lim l(Tr#4) ^8 - I6O5 + 4032ff2 + (2.5.3) 

N—*oo iv 
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whose first two terms stera from tlie planar graphs 



(2.5.4) 

^^^^X^ -325 Q -256g r ; 1285 S : -I6O5 

Planar diagrams with 2-fcrmionic legs are counted with the help of the observed 
Ward identities of eq. (2.1.20). For our quartic potential we have 

(Tr<P'<f') = i (Tr#2) + | (Tr^*). (2.5.5) 

Hence with the results of eqs. (2.5.1) and (2.5.3) the sum of the planar amplitudes 
containing 2-fermionic legs reads 

lim ^ (Tr^^) = 1 - 4g + 6Qg^ + (2.5.6) 

N^oo N 

and may be seen to originate from the diagrams 

1 E : 1 



-Ag 



(2.5.7) 



I652 



6452 _32g2 S : eOy^ 



For 



confirming our calculations. 

Of course our results allow the counting of non-planar diagrams as wel 
example the toroidal diagrams with 2-bosonic legs are contained in G\{g)^ i. 

G\{g)^ lim iV(Tr = _8g + 560g2 + _ (2.5.8) 
whose first term stems from the toroidal diagrams 

-I65 ^ 8g S:-?,g (2.5.9) 

This concludes our analysis of the O-matrix model. 
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6. Conclusions 

In this chapter we have studied the O-matrix model, defined by a matrix valued 
superfield in (0,2)-dimensional superspace. After integrating out the auxiliary and 
fermionic degrees of freedom its partition function was shown to be identical to 
the partition function of a Hermitian matrix model with pure Gaussian potential. 
These two models are related to each other through a Nicolai-map. The analysis 
of the Feynman rules of the O-matrix model revealed that its dual graphs may 
be interpreted as discretized random surfaces decorated with sign producing loops. 
Subsequently the one point-functions of the O-matrix model with quartic potential 
were calculated for all genera, as well as the planar two- and three-point functions. 
Despite of the triviality of the partition function, these correlators are non-trivial 
and display critical behaviour. The scaling limit of these correlators reproduced the 
string susceptibility exponent of c = — 2 Liouville theory. Finally we showed how 
these results represent the solution to the combinatorial problem of supersymmetric 
graph counting. 



Chapter III: 

The Supereigenvalue Model 



Following the successful application of Hermitian matrix models to 2d gravity 
and lower dimensional bosonic strings, it was natural to ask how these methods 
could be generalized to the supersymmetric case. Unfortunately the description 
of discretized 2d random surfaces by the Hermitian matrix model to date has no 
analogue in terms of some supersymmetric matrix model describing a discretization 
of super-Riemann surfaces. In a sense that might not be too astounding as there 
exists no geometrical picture of superspace. 

Nevertheless the supereigenvalue model proposed by Alvarez-Gaume, Itoyama, 
Manes and Zadra [30] appears to have all the virtues of a discrete approach to 
2d quantum supergravity. Guided by the prominent role the Virasoro constraints 
(cf. eq. (1.3.2)) played for the Hermitian matrix model, the authors constructed 
a partition function built out of N Grassmann even and odd variables (the "su- 
pereigenvalues" ) obeying a set of super- Virasoro constraints. In order to formulate 
these constraints it is necessary to augment the bosonic coupling constants gk by 
a set of anticommuting couplings usually denoted by S,k+i/2- As it is unknown 
whether there exists a matrix-based formulation of this model, we do not have a 
geometric interpretation of it at hand. Nevertheless many of the well known features 
of the Hermitian matrix model, such as the genus expansion, the loop equations, 
the double scaling limit, the moment description and the loop insertion operators, 
find their supersymmetric counterparts in the supereigenvalue model. From this 
point of view the supereigenvalue model appears as the natural supersymmetric 
generalization of the Hermitian one matrix model. 

The supereigenvalue model is solvable nonperturbatively in coupling constants 
but perturbatively in its genus expansion. Again the problem of solving the model 
may be reformulated in a set of superloop equations obeyed by the superloop cor- 
relators. Away from the double scaling limit these equations were first solved for 
general potentials in the planar limit in ref. [31]. We were able to develop a complete 
iterative solution of these equations for all genera and general potentials in ref. [5] . 
The key point in this scheme is the change of variables from coupling constants to 
the so-called moments of the bosonic and fermionic potentials, thus generalizing the 
approach of ref. [8] for the Hermitian matrix model. An alternative approach was 
pursued by the authors of refs. [32, 33] who managed to directly integrate out the 



22 



□ 

Grassmann odd variables on the level of the partition function. No generalization of 
a solution based on orthogonal polynomials is known. The supereigenvalue model 
also displays a connection to supersymmetric integrable hierarchies [34] . 

In order to make contact to continuum physics the supereigenvalue model has to 
be studied in its double scaling limit. This was done in refs. [SO? 35, 31] for pla- 
nar (and partially for higher genus) topologies. These studies revealed that in its 
continuum limit the supereigenvalue model describes the coupling of minimal su- 
perconformal field theories (of type (2, 4m)) to 2d supergravity. Moreover Abdalla 
and Zadra [36] were able to develop a dictionary between N — 1 super-Liouville 
amplitudes and supereigenvalue correlators. It turns out that the moment descrip- 
tion is very useful for the identification of the critical points in the space of coupling 
constants. In ref. [6] we applied this method to the scaling limit, and presented an 
improved iterative scheme to directly compute only the higher genus terms relevant 
in the double scaling limit. 

In the following we derive the supereigenvalue model directly from the super- 
Virasoro constraints in a superconformal field theory formulation of the problem. 
After some short remarks on the general structure of the free energy, we intro- 
duce superloop correlators and the superloop insertion oper ^tci| rs. The superloop 
equations are derived and solved through the iterative procedure discussed above. 
We then turn to the double scaling limit, identify the scaling of moments and ba- 
sis functions and present our improved iterative scheme. Finally some remQs on 
the identification of the double scaled supereigenvalue model with N — 1 | sij per- 
Liouville theory are made. 



□ 



1. Super— Virasoro Generators 



As we have discussed, we wish to regard the super- Virasoro constraints feHhe fun- 
damental properifti-of a supersymmetric generalization of the Hermitian one matrix 
) generators of the Hermitian matrix model in eq. (1.3.3) are dif- 



model. The Viraf or 



ferential operatorsm the coupling constants gk ■ In order to write down an analogous 
set of super- Virasoro generators it will be necessary to introduce anticommuting 
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coupling constants ^k+1/2 as the "superpartners" of the bosonic coupUngs gk- The 
super- Virasoro generators are then given by [30] : 

1 " 



k>0 



k>0 



k=0 



1 " ?l + 1 



/c>0 



fc>0 



n— 1 



2 7V2 



n > -1. 



(3.1.1) 



fc=0 



Note that in the above we define J^k^n ■ • • = 0. The super- Virasoro generators 
satisfy 

[ , ] = ^ Tl) Ln+m 

[^n , Q7n+l/2 ] = in/2 - 1/2 - to) tJ„+„i+i/2 
{5n+l/2 : Sm+1/2} = 2 (3.1.2) 

which constitutes a closed subalgebra of the = 1 superconformal algebra in the 
Neveu-Schwarz sector (remember that n,m> —1). 

2. Construction of the Model 

The supereigenvalue model will now be constructed from the requirement that its 
partition function Zs[gk,£,k+i/2\ is annihilated by the {§f^i/2T ^n] n > How- 
ever, equation (3.1.2) shows that it is sufficient to impost 

Gn+i/2Zs = 0, n>-l, (3.2.1) 

since the constraints iL„ Zs = for n > — 1 come out of the consistency conditions 
of eq. (3.2.1) using the algebra (3.1.2). 

The method to construct -Zs[ff/c: Cfc+i/2] is based on a comment made in ref. [37] 
to generalize the conformal field theory formulation of matrix models discussed in 
chapter I. We seek for a representation of Zg[gk, 'Cfc+1/2] in the form of a correlator 
in a superconformal field theory. Let us consider a free, holomorphic superfield 
X(z, 9) — 4>{z) + 9ijj{z) with the mode expansions 



(z) Inz - ^ ^ z 



i^{z)^^$k+l/2 

feez 



fc/O 



(3.2.2) 



3.2. Construction of the Model 
whose modes obey the (anti)-commutation relations 

[Jn,Jrn]=nSn+m,0, [?,?] = !, 

(3.2.3) 

{'0n+l/2 , 1?m-l/2} = Sn+m,0, [Jn , V'm ] = [q , V'm ] = [p , i'rn] = 0. 

Define the vacuum states 

Jk |0) = (27V| J_fc = k>0 
p\0)=0 

V^„+i/2|0)=0 (27V| ^_„_i/2 = n>0, (3.2.4) 

□ 

where |2A'') = exp[2iVg] |0) and q is anti-Hermitian. The super-energy-momcntum 
tensor of this theory is given by 

T(z, 61) ^ T^{z) + 2 0T'^{z) = ^d(l) + : {d(j)d(j) + ijdip) (3.2.5) 

whose fermionic part has the mode expansion 

T^(z) = ^rf^i/2Z— 2 (3.2.6) 

Tn+l/2 = ^ij-k'!pk+n+l/2+'!p-k+l/2Jk+n]+ ^ V'a+l/2^fc, n > 0. 



fc>0 



a-\-b—n 



Similar to the bosonic case we introduce a Hamiltonian built from Grassmann even 
and odd coupling constants 



Higk,^k+i/2)^^(^gkJk +6+1, 



k>0 



: 'ipk+i/i 



Co 



^(f(^)5(/.(z) + vI/(z)V-(z)), (3.2.7) 



where V{z) — ^k>o9k *(-^) 
shows that 

^n+i/2 {2N\ exp[iJ(gfc, 6+1/2)] 



'l2k>o^k+i/2 ■ Using these definitions one 



(2iV| exp[i/(5fc, 6+1/2) ]Tf+i/2 



(3.2.8) 
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III: The Supereigenvalue Model 
Hence any operator satisfying ^ 



-1/2 



O] = 0, 



n > -1, 



will give us a correlator obeying the super-Virasoro constraints 

5n+l/2 im exp[i7(5fe,Cfe+l/2) ] 0, 



(3.2.9) 



(3.2.10) 



as T^^-y^^ |0) =0. Thus every nonvanishing correlator of this form may be used for 
the definition of the partition function 



Zs = (2iV|exp[i7(5fe,efc+i/2)]O|0). 



(3.2.11) 



The problem of finding the operators obeying eq. (3.2.9) is solved in superconfor- 
mal field theory. It turns out that any function of the super screening charges Q± 
will do 



<j> dz j de : exp[±X(z,6i)] 



(3.2.12) 



where X(z, 6) again denotes the holomorphic superfield. We choose to get a 
nonvanishing correlator in eq. (3.2.11) ^ . Note that we take an even 27V in order to 
have a Grassmann even partition function Zs- The model is then defined by 

2's[gk,S,k+i/2 ] = 

2N 

{2N\ :exp[i7(5fe,Cfc+i/2)] Hf ^-^^ d0, : exp[(j){zi) + e,ij{z,)] : |0). (3.2.13) 

This correlator may be evaluated by using the operator product expansions 
(j>{z) 4>{uj) ^ ln(z — uj) and ipiz) ip(oj) ~ (z — ui)^^ ■ The result is 



-2s [.g*;, 4+1/2] 



2N 2N 

Wj dz, / exp[^(l/(zO + ^^*(^0)] n( 



,). (3.2.14) 



Note the resemblance of this model to the eigenvalue description of the Hermitian 
matrix model of eqs. (1.1-4) and (1.3.16). This is the reason for calling Zs the 



^ We again extend the definition of T-^, , ,„ 

n+1/2 

second eq. of (3.2.6) 

^ More generally one could take Q 
have not been analyzed so far. 



to 71 = —1 by dropping the second term in the 
i2JV+A/ qM ^ supereigenvalue models constructed in this fashion 



3.3. Structure of the Free Energy 

"supereigenvalue model" despite of the fact that the Zi and 9i are no eigenvalues of 
any matrix It was first constructed in ref. [30] using different methods. 



The derivation of eq. (3.2.14) goes along the same 1: net as the bosonic case of eq. (1.3.16). Using 
: : : : = : e^+B we have 



2N „ „ „ 



2JV 



Co 



i = 2 
2JV 



:exp[fli(z,) + 6»ii/.(zi)] : |0> • - |__| 

j> dz, / de,^ {2N\ : exp[ (L [ V(z) d<i>(z) + m{z) iP{z) ] + ^ ^(zj) + 9^ V(2j) ] : 
Jl :exp[fli(z,) + f»^(z,)] : |0> e'"'"^"""'"^!^ exp[^ y(^,) + 6^ vl>(^ ■)] 



i=3 



j = l 



2JV „ „ „ 2JV 

(jlj dz, j de^{2N\:cM<f [V(z)d<t>(z) + 'i{z)^l)(z)]+^<t>{z,) + e,^(z,)]:\fi) 

2JV 

■ - z, - Q,e,) cxp[^ y(^j) + ^{z^)\ 

i<j j=l 

2JV „ „ 2JV 

(n f / '^^') n^"^' ~ ~ ^''^^^ exp[^ V{z,) + e, ^{zj)] ^(2Afh_cxp [ 2Afg] : |0)^ 



and we recover eq. 3.^ .14) 



3. Structure of the Free Energy 

After deforming the contours of integration in eq. (3.2.14) to the real line we recover 
the supereigenvalue model proposed by Alvarez-Gaume, Itoyama, Mafiis and Zadra 
[30]. After a slight change of notation the partition function is given ay 



''Thc only attempt in literature to formulate eq. (3.3.1) in terms of matrices may be found in 
ref. [38], where an explicit solution to this problem was found only for low N. 
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Zs[9k,£.k+i/2;N] = 

2N 



3.3. Structure of the Free Energy 
where Zsi'^gk', N) is the partition function of the purely bosonic N x N Hermitian 



2N „oo n 2Af 

[](/ dX, / l[{X,-X,-e,9,)exp[-2NJ2[y\^^)-S^'fi^^))], (3.3.1) 



1<J 



matrix model of eq. (1.1.1). We pr 

□ 



built from a set of 2N bosonic and fermionic variables and 9i respectively. The 
bosonic and fermionic potentials read 



V{Xi}=Y,9kX,\ and *(A,) = ^ efc+i/2 A,^ 



(3.3.2) 



k=0 



k=0 



with the Grassmann even and odd coupling constants gk and ^fe+1/2 respectively. 
The free energy of the supereigenvalue model is naturally defined by 



25(<?.,6+i/2;27V) = e(2^)^^(^-«^+v.;2A'). 



Before developing an iterative procedure to solve the supereigenvalue model, let 
us discuss the fermionic structure of the free energy. Consider the supereigenvalue 
model of eq. (3.3.1) in the absence of fermionic couplings, i.e. S,k+i/2 = 0. The 
integration over the fermionic variables 6i may then be performed ^ giving rise to 
the effective bosonic eigenvalue model 

Zs[gk,0;2N] = 



27V 



2N 



n( / '^^^ ) n(^' - )Pfaff (A-.i) exp ~2Nj2viXi) , (3.3.4) 



where Pfaff (A^^^) is the Pfaffian of the antisymmetric matrix 



A, - A, ' 



Pfafr(A) ^VdetT^ 



1 



2^iV! 



(3.3.5) 



There exists an interesting identity between the effective bosonic model of eq. (3.3.4) 
and the Hermitian matrix model 



Zs(.gfe, 0; 2N) = c{N) ZB{2gu; N) 



(3.3.6) 



*Simply write n.<,(^' " " ^»«.) = U^<,(>^^ " ^ME.<, xf=iq] 



The statement wc wish to prove is 

/2JV 
(j[ <iM(A«)) PiaSiXr^"-) A(Ai, . . . , Aziv) 



)vt this identity in the following. 

□ 



(N) / rj^dMA,)) A2(Al,...,A^r)A2(AJv+l,...,A2Jv), 



(3.3.7) 



where /i(Ai) is a measure factor depending only on the eigenvalue A^, A(Ai, . . . , Xm) denotes the 
van der Monde determinant and c{N) some irrelevant TV dependent factor. 



We prove eq. (3.3.7) by induction in A'^. First note that due 
A(Ai, . . . , A2jv) under the exchange of two eigenvalues we have 



/2JV 
(j[ dfj.(\,)^ Pfaff(A-:i) A(Ai, . . . , X2n) 



he antisymmetry of 



zi\ I 

ill' 



dii{Xi) A(Ai,...,A2]v) 



Ai — A2 A3 — A4 A 



2JV-1 — ^2N 



(3.3.8) 



up to irrelevant factors depending on TV. Moreover from A(A) = (— )'^(^ l)/^ det(Ai.' ^) we have 
the identity 



]^dM(A,)) A2(Ai,...,Ajv) = / (^J|rfMA>)) A(Ai,...,Ajv) J|a,'-i. 



(3.3.9) 



1=1 1=1 
Let us now show that eq. (3.3.7) is true for TV = 2 
4 



Za = j dM(A«)) (Ai - A3) (A2 - A3) (Ai - A4) (A2 - A4) 

i = l 

= / (n "^^(^'0 '^"^ (1 A3 ) ( 1 A4 ) ~ ~ 



i = l 
4 



det ( I \ 1 det 
A3 A3 



the assumption (3.3 



the second term in the last expression vanishes. Applying eq. (3.3.9) to the first term then proves 



7) for TV = 2. 



Now by using the induction hypothesis the left hand side of eq. (3.3.7) for (TV + 1) becomes 
2JV+2 2JV 



Y\ dix{Xi)^ A2(Ai, . . . , A^r) A2(A^r+l, . . . , A2iv) J|(A, ■ 

_ 1 „■ 1 
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2JV+2 



Jl dA<(AOj A{Ai,...,Ajv) J|a/-i A{Ajv+i,...,A2jv) A, 



2JV 



□ 

• ]^(Ai — A2]V+l) (Ai — A2JV+2) 

i = l 

„ 2JV+2 JV 
- ^M^o) A(Ai,...,AAr,A2iv+i) J|a,'-i A{Ajv 

i—1 i—l 
2N N 

i+JV - A2JV+1) (Xi - A2JV+2)- 



j=JV+l 



, A2]V, A2JV+2) 



(3.3.10) 



i=JV+l 



Note that 



A(Ai, . . . , A^r, A2]v+i) Ai' 



: det 



/ 1 Ai Ai2 

' A2 A22 Aa^ 



\ 1 A2JV+1 A2]V+1^ ... A2]V+1^^^ / 



as well as an analogue expression for the second van der Monde determinant in eq. (3.3.10). One 
may then convince oneself that in the product 



J^(Ai+JV — X2M+1) (Xi — X2N+2) — 

i = l 

N 

EE E (-)'^-+'wi^-'=A,, 



• Ai, \2N+2^ ' Xji ■ • • Aj^, 



e{l,...,N) G(Af + l,...,2Af) 



only the following terms contribute to Z2JV+2 of eq. (3.3.10) 

JV 

A2JV+l"'^ A2JV+2"'^ + ^ ^ A2JV+1'" A2]V+2' Ar A^+l . 



. A]V Ajv+r Ajv+r + 1 . • • A2JV 



A2 



Hence after some relabeling of indices and via eq. (3.3.9) we arrive at 

/2JV+2 JV 
[| dM(A,)) A(Ai,...,Aiv,A2iv+i) H^*'"'^ 
i—l i=l 
2N 

A(A]v+l, . . . , A2iV, A2]V+2) Xi"^^ ^ X2N+2 

. iJV+2 

( Jl dM(A,)) A2(Ai,. 



2JV+1 



i=JV+l 



Ajv,A2jv+i)A (Ajv+1, . . . , A2]v, A2JV+2), 



□ 

3.4. Superloop Equations 

and we have proven eq. (3.3.7). ' ' 

Eq. (3.3.6) was first proposed in ref. [32]. It implies that up to irrelevant additive 
constants the part of the free energy independent of the ferniionic couplings obeys 

F^s^[gk,^k+i/2 - 0;2iV] = 2Fs[fffe;iV], (3.3.11) 

a relation which we will recover in the explicit solution of the supereigenvalue model. 

Including the fermionic couplings Cfe+1/2 in 6q. (3.3.1) and performing the 
fermionic integrations to obtain the complete effective bosonic eigenvalue model 
is a more complicated task. Nevertheless McArthur [33] succeeded in doing so. 
The outcome important to us is the remarkable result that the free energy contains 
contributions only up to second order in the fermionic couplings, i.e. 



Fs[g^k+i/2] = i^f [5fc] + E^'«+i/2?'+i/2^M [f'c], (3.3.12) 



k+l/2\ 

k,l 



confirming a conjecture made in ref. [32]. This observation will be the key to the 
solution of the superloop equation which wc will describe in the following section. 



4. Superloop Equations 

Our solution of the supereigenvalue model is based on an integral form of the super 
loop equations, which are the analogue 
trix model discussed in chapter I. Gene: 



the loop enuations of the Hermitian ma- 
zing the kp Droach of Ambj0rn, Chekhov, 
Kristjansen and Makeenko [8] for the Hermitian one matrix model we develop an 
iterative procedure which allows us to calculate the genus g contribution to the 
(n|77i)-superloop correlators for (in principle) any g and any (n|m) and (in practice) 
for any potential. The possibility of going to arbitrarily high genus is provided 
by the superloop equations, whereas the possibility of obtaining arbitrary (n|m)- 
superloop results is due to the superloop insertion operators introduced below. A 
change of variables from the coupling constants to moments allows us to explicitly 
present results for arbitrary potentials. The remainder of this chapter is essentially 
an enlarged version of the authors papers [5] on the general and [6] on the double 
scaled solution of the supereigenvalue model. 

Superloop Insertion Operators 

For notational simplicity let us from now on write the supereigenvalue model with 
N taken to be even as 
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f{f[dX,d9,) [](A, - X, - e.Oj) exp(-7V^ [y(AO - f?.*(A,)]), (3.4.1) 



N N 



i—1 i<j i—1 



with the potentials V{Xi) = J2T=o dkXi' and *(Ai) = X^fclo ^^+1/2^4'' of eq. (3.3.2). 
Expectation values are defined in the usual way by 

(0(X,,6,)) = 

„ AT N 

^ {\{dx,de,) ^{x,,ei)0{x,,ej)e^^{~NY,[v{Xi)-e,^{Xi)]), (3.4.2) 
^ i=i i=i I I 

where we write A(Ai, 9i) — Y\i^j{Xi — Xj — OiOj) for the measure. We4ntroduce the 
one-superloop correlators 

W{p\)^nIJ2-^) and W{\p) = Nlj2-^y^ (3-4.3) 

i i 

which act as generating functionals for the one-point correlators {'Ylii^i') ^^^d 
( X^i ^i'^i'^ ) upon expansion in p. This easily generalizes to higher-point correlators 
with the (n|m)-superloop correlator 

W{pi, . . . ,Pn I 91, . . ■,qm) = 

^„+,n/y^^^_ v-^^y-^^ V- 1 \ (344) 

\^Pi- A,, ^ p„ - A,„ gi - Xj, ^ qra - Aj„ / 

Quite analogously to the bosonic case these correlators may be obtained from the 
partition function Zs by application of the superloop insertion operators 5/5V{p) 
and 5/5'^ij)): 



W{pi, . . . ,Pn \ qi, ■ ■ ■ ,qm.) = 

IS S S 5 



ZsS^ipi) 5^{pn) 5V{qi) 5V{qm) 

where 

5 ^ I d , 8 



Zs. (3.4.5) 



El c; 6 1 



5V{p) f^^p'^+^dgk 5m{p) f^^p^+^ d^k+i/ 



3.4. Superloop Equations 



In particular equation (3.4.3) can now be written as W{p \ ) = S lnZs/(5^(p) and 
W{\p) = 6\nZs/SV{p). 

However, it is convenient to work with the connected part of the {n\m)- 
superloop correlators, denoted by W. They may be obtained from the free energy 
F = N^^ In Zs through 



W{pi,...,p„ I gi,...,g™) : 
5 5 6 



□ 

F= (3.4.7) 



S-^{p,) ■ ■ ■ S^ip„) SV{q,) ■ ■ ■ SV{qm) □ 

\ ^ Pi - Ail ' ^ P« - Ai„ ^ 9i - Aji "' ^ qm~ Aj„ 

in complete analogy to the bosonic case of eq. (1.1.10). Note that n < 2 due to the 
maximally quadratic dependence of F on the mentioned in section 3. 

With the normalizations chosen above, one assumes that these correlators enjoy 
the genus expansion 

00 ^ 

W{px, . . . ,p„ I qi, . . . ,g„0 = ^ M/'s(pi, . . . ,p„ I gi, . . .,*«). (3.4.: 

9=0 

Similarly one has the genus expansion 



for the free energy. In contrast to the Hermitian matrix model there is no geo- 
metrical argument based on Feynman diagrams and the Euler relation available to 
justify the genus expansion. This is of course due to a missing formulation of the 
supereigenvalue model in terms of some generalized matrix model. However, the 
genus expansion is motivated from the structure of the superloop equations, which 
we derive in the following subsection. 

Superloop Equations 

The superloop equations of our model are two Schwinger-Dyson equations, which 
are derived through a shift in integration variables Xi and Qi. They were first stated 
in [30, 31], and we present them in an integral form for the loop correlators W{j) \ ) 
and W{\p). 

The Grassmann-odd superloop equation reads 
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r duo V'{uj) , , , r du ^{uj) . . 
(h — W(uj \ ) + i — W( uj) 



1 



W{p\)W{\p) + W{p I p) 



(3.4.10) 



and its counterpart, the Grassmann-even superloop equation, takes the form 
c2ti p-oj ^ ' ' Jc 2tt p-oj ^ ' ^ 2 dp Jc 2tt p - uj 



\ \w{\pf - W{p\)W\p\) + -^(w{\p,p) - (3.4.11) 

In the derivation we have assumed that the loop correlators have one-cut structure, 
i.e. in the limit N ^ co we assume that the eigenvalues are contained in a finite 
interval \x.y\. Moreover C is a curve around the cut. 



In order to derive the first of the two superloop equations for the model of eq. (3.4.1) consider 
the shift in integration variables 



\i — ^ + Qi - 



and 



where e is an odd constant. Under these we find that 



9^ + 



P- 



and the measure transforms as 



The vanishing of the terms proportional to e then gives us the Schwinger-Dyson equation 

1 



: 0. 



(3.4.12) 



Note that ( di(p — \i) ^ .(p — Xj) ^) = N ^ W{p \ p) with the definitions of the previous 
subsection. In order to transform eq. (3.4.12) into an integral equation we define the bosonic and 
fermionic density operators 

pW = ^^{K^-K)) and r(A) = i ^(e,5(A-A0>. 



With these the first sum in eq. (3.4.12) may be written as 



dXi r(A) 



p-X 



p- A 



/ dX r(A) 





1 


V'iu) 


2tt lu 


- A 


p — U) 



+ N'^ dX p{X) 



3.4. Superloop Equations 

du) 1 •^{ui) 



2-K u) — X p — u) 



(3.4.13) 



where we assume that the real eigenvalues A^ are contained within a finite interval x < A; <y 
in the limit N ^ oo. Moreover C is a curve around the cut [x, y\ and p lies outside this curve. Per- 
forming the A integrals in eq. (3.4.13) gives us the one— superloop correlators. The full Schwinger- 



Dyson equation (3.4.12) may tl 



du) V'(i^) 



2-K p ~ iO 



be expressed in the integral form 



^ti^L I) + (b — 



du) ^{(jj) 



2n p ~ ui 



Wi\uj) = N-^W{p\p). 



Rewriting this in terms of the connected superloop correlators W of eq. (3.4.7) yields eq. (3.4.10). 

The derivation of the second superloop equation goes along the same lines by performing the 
shift 

, 1 ef. 



A, 



Ai + ■ 



and 9i 



p-Xi " 2(p-A,)2 
with e even and infinitesimal. Similar steps as the ones discussed above then lead us to the second 
superloop equation 



C 



dw V'juj) Wi \ w) + ^'{uj)Wiui I) 
2tt p — uj 



1 d 

2 dp 



1 1 



W(\p,p) 



c 

W{p,q\) 



dio ^{u))W{u) I) 
2-K p — Ul 



lid 

2 " " 2^ dq 

which, after rephrasing in connected quantities, gives eq. (3.4.11) 



Note the similarity to the loop equation for the hcrmitian matrix model (1.3.21). 
The eqs. (3.4.10) and (3.4.11) are equivalent to the super- Virasoro constraints. 
The constraints follow from the superloop equations by expanding them in p and 
evaluating the resulting contour integrals. The resulting equations at each power 
of p then correspond to Gn+i/2 — and Zs — for n > — 1. 

The key to the solution of these complicated equations order by order in N^"^ 
is the observation stated in eq. (3.3.12) that the free energy F depends at most 
quadratically on fermionic coupling constants [32, 33]. Vwi'^q- (3.4.7) this directly 



translates to the one- loop correlators, which we from nov on write as 

W{p\)^v{p) 
W{\p)^u{p) + u{p). 



(3.4.14) 



Here v{p) is of order one in fermionic couplings, whereas u{p) is taken to be of 
order zero and u(p) of order two in the fermionic coupling constants £,k+i/2- This 
observation allows us to split up the two superloop equations (3.4.10) and (3.4.11) 
into a set of four equations, sorted by their order in the (,k+i/2^^- Doing this we 
obtain 



Order 0: 
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dio V'ico) 

U(UJ) — 

ZTT p ^ UJ 

1 . ^2 11 S 



lid S 



Order 1: 



(3.4.15) 
□ 



du V'iu;) 
' t; ^\^) 

C P — U! 



duj *(a;) , , 

u{uj) 

(J 2tt p — iv 



v{p) u{p) 



1 



7V2 SV{p) 



v{p) 



(3.4.16) 



Order 2: 



dio V'iuj) 
2ti p ^ bJ 



u{lo) + 



duj *'(a)) 



u{p) u{p) 



c 2tt p 
1 



1 d 

2 dp 



duj *(u;) 
27r p — 



2 Tp "(P) 



1 1 



2 iV2 5V(j)) 



(3.4.17) 



Order 3: 



duj '^/{uj) 

u{uj) = v{p)u{p). 

Q Ztt p — uj 



(3.4.18) 



It is the remarkable form of these four equations which aUows us to develop an it- 
erative procedure to determine Ug{p), Vg{p),Ug{p) and Fg genus by genus. Plugging 
the genus expansions into these equations lets them decouple partially, in the sense 
that the equation of order at genus g only involves Ug and lower genera contribu- 
tions. The order 1 equation then only contains Vg, Ug and lower genera results and 
so on. The first thing to do, however, is to find the solution for <? = 0. 



5. The Planar Solution 

In the following the planar solution for the superloop correlators is given for a 
general potential. It was first obtained in [31]. We present it in a very compact 
integral form augmented by the use of new variables to characterize the potentials, 
the moments. 
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Solution for uo{p) and vq 



P) 



In the limit iV — > oo the order equation (3.4.15) becomes 



duj V'{lu) 
C 2t: p — uj 



ua(uj) 



■uo{pf 



□ 



(3.5.1) 



This equation is well known, as up to a factor of 1/2 it is nothing but the planar 
loop equation of the hermitian matrix model (1.4.2). With the above assumptions 
on the one-cut structure and by demanding that u(p) behaves as 1/p for p — > oo 
one finds [16] 



uo{p) 



du V'{lu) 
U 2t: p ~ lu 



{p-x){p-y) 



1/2 



(3.5.2) 



(w - x){u) - y) 

as we verified in chapter I section 4. The endpoints x and y of the cut on the real 
axis are determined by the following requirements: 

duj V'{uj) r duj ujV'{uj) 







c 271- ^{uj ~ x){u} -y)' 



1 



c 2vr ^ [lo - x){uj - y)' 



(3.5.3) 



deduced from our knowledge that M^(|p) = l/p + 0{p^^). 

The order 1 equation (3.4.16) in the N ~t oo limit determining the odd loop 
correlator vq{j)) reads 



duj V'{uj) , , r dcj *(w) , , , , , , 

vo{uj) + f Mo(w) = vo{p)uo(p). 

Q In p — LO Jfj In p ~ UJ 



It is solved by 



□ 



duj *(w) 
P 27r p — 



{uj - x){uj - y) ] 
{p-x){p-y) 



^/{p-x){p-y) 



(3.5.4) 



(3.5.5) 



Here % is a Grassmann odd constant not determined by eq. (3.5.4), in fact 
X = N^^ ( J2i ) ill planar limit. It will be determined in the analysis of the 
two remaining equations (3.4.17) and (3.4.18). 

One verifies the above solution by direct computation. Using eqs. (3.5.2) and (3.5.5) the right 
hand side of eq. (3.5.4) reads 



uoip) vo{p) ■- 



du) 
2^ 



2tt (p - ui){p - z) 
X 



{z - x)(z ~ y) 



Ci JC2 

2n p — UJ [(uj ~ x)(ll} ~ y)]^/'^ ' 



(w - x){uj - y) 



1/2 



(3.5.6) 
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3.5. The Planar Solution 



and the left hand side becomes 
duj f dz y'(a;)<^(z) 



Ci Tc2 ^'^ ~ '^^('^ ~ 



(z - x){z - y) 



-I 1/2 



(a; - x){u; - y) 
duj f dz ■^{ui)V'(z) 



{iv - x){uj - y) 



1/2 



The main motivation for introducing these new variables is that, for each term in 
the genus expansion of the free energy and the correlators, the dependence on an 
infinite number of coupling constants arranges itself nicely into a function of a finite 
number of moments. 

We further introduce the basis functions x'""'\p) ^^^d ^''^"^(p) recursively 



(z - x)(z - y) 

Now in the last term pull the contour integral C2 over the curve Ci. One can show that the 
contribution from the extra pole vanishes. After renaming ui ^ z and combining the first and 
third terms one gets eq. (3.5.6). Thus eq. (3.5.4) is verified. 

Moments and Basis Functions 

Let us now define new variables characterizing the potentials V{p) and ^{p). Instead 
of the couplings gu we introduce the bosonic moments Mk and Jk defined by [8] 



ri-fc+l 



whi 



dw V'{uj) 



C 27r {oj-xY [(w-a;)(w-y)]i/2^ 
1 



k > 1 
k>l, 



and the couplings £,k+i/2 a-re replaced by the fcrmionic moments 



(3.5.7) 



(3.5. 



following ref. [8]. 

It is easy to show that for the linear operator V defined by 

V'ofip)^fp^fiu)-uoip)fip) 
Jc 271 p-w 

and appearing in the superloop equations we have 

1 



{p - xY 



n > 1, 



These moments depend on the coupling constants, both explicitly and thfou|gh x 
and y: 



Mfc = (/c + l)gfc+i+ (fc + i)x + izj {k + 2)gk+2 + 



Jk^{k + l)gk+ ^x + {k+^)y (fc + 2)5fc+2 + ... 



.^fe = ( 1 - (5i,fc ) ^fc_3/2 + {k-lr)x-l-y + 



Afc = ( 1 - ) ^fc_3/2 + {k-^)y-^x + 



□ 



(3.5.9) 



V' o vl,(")(p) = _^_, n>l. 

Moreover, (t>l^\p) = 4'y'\p) = 't'^^Hp) lies in the kernel of V. 

Solution for uq and x 
Next consider the order 2 equation (3.4.17) at genus 

ex. ^ 1 d f du! *(w) , , f duj *'(w) , , 
^ dp J c 27r p — a; J(j Zn p — uj 

-^Mp)-^Mp)- 



(3.5.10) 



(3.5.11) 



(3.5.12) 



(3.5.13) 



(3.5.14) 
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Plugging eq. (3.5.5) into the right hand side of this equation yields after a somewhat 
lengthy calculation 



2 {x ~ y) p — X 



1 A2 (Ai - x) 1 

2 (a;-y) p-y' 



(3.5.15) 



(3.5.16) 



With eq. (3.5.13) this immediately tells us that 

There can be no contributions proportional to the zero mode (p), as we know 
that u{p) behaves as 0{p~'^) for p — > oo. 

Finally we determine the odd constant x- This is done by employing the order 3 
equation (3.4.18) for g = 0, i.e. 



uo{Lo) - vo{p)uo[p) =0. 

(J 2tt p — oj 

After insertion of eqs. (3.5.5) and (3.5.16) one can show that 

. ^ ^ , , 1 S2 (Si - x) (Ai - x) 1 A2 (Ai - x) (Si - x) 
"^P' 2 Ah{x-ynp-y) 2 J, {x - y^ip ~ x) ' 

where we have defined the linear operator V by 

duo *(w) 



¥ o f{p) 



f{uj) -vo{p)f{p), 



(3.5.17) 



(3.5.18) 



(3.5.19) 



j(j2TTp—UJ 

in accordance to V. The result (3.5.18) lets us finally read off the coefficient x as 

(3.5.20) 



1 



X = - (Si + Ai). 

Putting it all together, we may now write down the complete genus solution for 
the one-superloop correlators W{ \ p) and W{p \ ): 



Wo{\p) 



4 Mi{x-y) 



{p-x){p-y) 
(w - x){uj - y) 

(P) 



-.1/2 



du V'ico) 
Q 2Tr p — uj 

1 S2 (Si - Ai) ^(1) ^ 1 A, (S^Ai, ^(1) 



4 JiC 



Wo{p\) 



did ^{lu) 

(J 271 p — LO 



{uj ~ x){uj - y) ] 
{p-x){p-y) 



(3.5.21) 



Ai 



2 [(p-x)(p-y)]i/2- 
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One can show that this solution is equivalent tEOhe less compact one obtained in 
ref. [31]. |— I 

We shall make use of the following rewriting of the purely boson ic jart of Wo( | p) 



Uo 



^ OO 

ip) = V'ip) - - [{p-x){p-y)]'^' {(P-^)"'' + ip-yy-' J,}, (3.5.22) 



9=1 



derived by deforming the contour integral in eq. (3.5.21) into on; surrounding the 
point p and the other enifir^ling infinity. To take the residue at innmty one rewrites 
(p — uj)~^ as 



1 



1 



1 



p-LU 2 {p- x) - {ui - x) 2 {p- y) - {uj - y)'' 



(3.5.23) 



and expands in (^^rrf j and {^zj^J I'espectively. Doing this for the fcrmionic Wo{p \ ) 
yields 



ip) = 2 iiP-^)(p-y)]^'^' E {(p-^r-'^g + {p-yy-'Ag}. (3.5.24) 



q=2 



It is important to realize that the bracketed terms in eq. (3.5.22) as well as in eq. 
(3.5.24) are actually identical. Here we see that the planar solution is special in the 
sense that it depends on the full set of moments. Interestingly enough this is not 
the case for higher genera. 



6. The Iterative Procedure 



lution of the superlo(» 



equations results in a certain representation 



Our iterative sc 

of the free cne-gy and the loop cotrelators in terms of the moments and basis 
functions defined in section 5. We wilLslhow that it suffices to know Ug{p) and Vg{p) 
only up to a zero mode in order to calculate Fg. We give explicit results for genus 
one. 



The Iteration for Ug and Vg 



The correlators Ug{p) and Vg{p) are determined by the order and order 1 equa- 



tfOI] 

pe 



s (3.4.15) and (3.4.16) after insertion of the genus expansions (3.4.i 
ators. We find 

1 



of these 



V' 



-,{P) = 



9-1 



Ug'(jp) Ug^g'ip) 
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1 5 



2 SV{p) 



1 d 



2 dq S-^ip) 



p=q 



(3.6.1) 



3.6. The Iterative Procedure 



for fc > 1, one finds 



and 



d 5 



5V{p) 



ua{p). 



V O Vg{p) ^ O Ug{p) + ^ Vg,{p)Ug^g,{p) 



SVip) 



'^g-^ip) (3.6.2) This enables us to write down It i(p), 



at genus g > 1. From the structure of these equations we directly deduce that Ug{p) 
and Vg{p) will be linear combinations of the basis functions x^^^\p) f^^d 'I'(")(p). 
By eq. (3.5.13) the coefficients of this linear combination may be read off the poles 
{p— x)~'' and {p — y)"^ of the right hand sides of eqs. (3.6.1) and (3.6.2) after a 
partial fraction decomposition. 

Let us demonstrate how this works for g = 1. According to eq. (3.6.1) for ui{p) 
we first calculate 5uq/5V{p). We then need to know the derivatives 5 x/5V{p) and 
5y/5V{p). They can be obtained from eq. (3.5.3) and read 



lip) = \x^'\p) + W'\p) - ix^'^ + i-M'\p) 



8 



8 



Ad- 



Ad 



(3.6.8) 



(3.6.9) 



5x _ 1 
5V{p) ^ Ml 



Using the relation 



5V{p) 



5v{p) Ji 



Note that up to the overall factor of two this is identical to the one-loop correlator 
of the hermitian matrix model [8], as it has to be due to eq. (3.3.11). 

Now we solve eq. (3.6.2) at g — 1 for vi (p) . It is important to realize that generally 
eq. (3.6.2) fixes Vg{p) only up to a zero mode contribution Kg (j)^'^'^ (p). This comes 
from the fact that, unlike for the bosonic u{p), we do not know the coefficient of the 
p~^ term for v(p). The zero mode coefficient Kg will be fixed later on by r equ iring 
v „ (p ) to be a total derivative of the free energy Fg 



(3.6.3) 



Ir order to calculate 5vo/5V{p) we make use of the relation 



' cqu i: 



1 



dp p — U! 



one finds ^ 
6 



SV{p) 



ut){p) = 



8{p-x)'^ 8(p-j/)2 Ad{p-x) Ad{p-y) 



(3.6.4^ 



(3.6.5) 



5V{p) 



+ 



1 ^ 



2 {-dY+i^^^ 

r=2 ^ ' 



{-dY 



(3.6.10) 



as well— as 8t^kl^yij>) obtained from the above by the replacements x 
Mk ^ Jk 
calculate 



<-> Kk and d —d. The derivatives 5 Mk/ 5V{p) and 5 Jk/5V{p) were 
a in ref. [8] 



where d — x ~ y. 

Next we determine S vo{q) /S'i>{p). Using the relation 



SMk 



5V{p) 



" -\ ip - x)-^-^'\p - y)-^/2 -{k + 1/2) 4'=+i)(p) 



J*(p) 



*(<z) - - 



1 



(3.6.6) 



1 



E 



1 



M,, 



and the result 



(jj — xY 

[{p-x){p-y)Y/^ 



-(fc + 1/2) 



M, 



fe+i 



Ml 



(3.6.11) 



(5*(p) 
(5*(p) 



3fcl 



{p - y Y 



(3.3.7) 



and 5 Jk/5V{p) is obtained by the usual replacements. Using these and the earher 
results one has 



^Of course this is nothing but 2 Wo(p,p) of the Hermitian matrix model of eq. (1.4.7). 



5vo 

5V{p) 



] 1 


\ (Si-Ai)i 


1 


(p — xY 


V AdJi \ 


{p - y Y 
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3.6. The Iterative Procedure 



AdMi 

A2 
4d2 Ji 

^2 



1 



" 4^2 Ml 
A2 



A2 
4d Ji 



1 



4d2Mi 4d2Ji 



{p - x) 
1 

{p~y)' 



ip - vY 



For the evaluation of the right hand side of eq. (3.6.2) at ge luf 



(3.6.12) 

g we also need to 



know how the operator ^ acts on the functions (jrx (p) and— ^1, (p), in terms of 
which Ug{p) is given. A straightforward calculation yields 



Tl+l 



1 



1 



(p-y) 



2(„d)"+2-fc 

(Si-Ai) 



E 

1=2 



(3.6.13) 



2(-d)"+i 

as well as the analogous expression for W o 0y"''(p) obtained from eq. (3.6.13) by 
the replacements x ^ y, Mk ^ Jk, Sfc ^ Ak and d — d. 

We now have collected all the ingredients needed to evaluate the right hand side 
of eq. (3.6.2). After a partial fraction decomposition we may read off the poles at 
X and y, and therefore obtain the coefficients of the linear combination in the basis 
functions. We arrived at the result for g — 1 with the aid of Maple, na mely 
3 

z;i(p) =^(b«X^*^(p) + e['^^>(^\p)) + AJi0W(p), 1-1(3.6.14) 

1=1 

where the coefhcients b[^^ and E'j'-' are given by 



(1) 



1 



1 Sa 



1 A, 



8 dMi 8 d'^Mi A d^Ji 
IM2S2 J_M2(Si-Ai) 
16 



8 dMi^ 
3 (Si-AQ 
' 16 d^ Ml 



d2 Ml 2 

_3_ (Si-Ai) 
16 



1 J2(Sl 

16 



Ai) 



d2 



d^Jl 



(2) _ 



1 S2 ^ 1 M2(Si-Ai) 
16 dMi^ 



3 (Si-Ai 



dMi 



16 d2 Ml 



^(3)_ 5 (Si-Ai) 
^ 16 dMi 



(3.6.15) 



and £^1'^ = ^(M J, S ^ A, d ^ -d). 

Yet Ki is still undetermined. To compute it and the remaining doubly fermionic 
part Ug [p) of the loop correlator Wg ( | p) one can employ the order 2 and order 3 
eqs. (3.4.17) and (3.4.18) at genus g. It is, however, much easier to construct the 
free energy Fg at this stage from our knowledge of Ug(p) and Vg{p). 

The Computation of Fg and Kg 

As mentioned earlier, the free energy of the supereigenvalue model depends at most 
quadratically on the fermionic coupling constants. In this subsection we present an 
algorithm which allows us to determine Fg and Kg as soon as the results for Ug[p) 
and Vg{p) (up to the zero mode coefficient Kg) are known. 

The purely bosonic part of the free energy Fg is just twice the free energy of the 
hermitian matrix model. By using the results of Ambj0rn et al. [8] one may then 
compute the bosonic part of Fg from Ug{p). Here one rewrites x^"-'(p) and ^("^(p) 
as derivative vith respect to V{p). One easily verifies the following relations for 



two lowest basis functions: 



X^'\p) = 
^(i)(p) 

X^'\p) 

*(2)(p) 



5v{py 



5y 
5V{py 

(5 2 1 
__(__lnj,__lnd). 



(3.6.16) 



Combining this with the result for ui{p) of eq. (3.6.9) yields the purely bosonic 
piece of Fi 

Fi^- = InAfi - 1 InJi - 1 Ind, (3.6.17) 

which as expected is just twice the result for the Hermitian matrix model (cf. [8]). 

For the higher basis functions the situation is not as simple. However, a rewriting 
of the basis functions allows one to identify Ug{p) as a total derivative. For x^"'\p) 
one uses the recursive form [8] 



Hi 



2n- 
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2 5 M„_ 



n-l 



2n - 1 6V{p) 



--Y.X^''\p)M^^u+i n>2, (3.6.18) 



k=2 



where 



should be rewritten as 

n-l 



2 5Af„_i 



4i)(p) = Mi 



(5 X 

5V{p) 2n - 1 SV{p) 
S X 



SVip) 
n> 2 



svip) ■ 



(3.6.19) 



The rewriting of ^'("^(p) foUows from the above by performing the usual replace- 
ments J ^ M and d ^ —d. 

The strategy for the part of Fg quadratic in fermionic couplings consists in rewrit- 
ing Vg{p) as a total derivative in the fermionic potential '^(p). We know that eq, 
(3.6.7) implies 

'^(Hi-Ai) 



(5*(p) 



and 



where 



J*(p) 



d(^\p) ^{p^x)-'^[{p^x){p-v)Y/^ k>l, 
<'Hp) = (p - 2/)-' [ {P - a;)(p - y) A: > 1. 



(3.6.20) 
(3.6.21) 

(3.6.22) 



Let us again specialize to 5 = 1. Using the above we can reexpress eq. (3.6.14) as 

3 



J*(p) 



Fi + Hi 



1 5 



d ^*(p) 



(Si - Ai) = E(^i'' ^^^^^P) + 



i=l 



Y^h^i^ ^^;\p) + 6('-)4'-)(p)) + 7i[4'Hp)-<)(p)], (3.6.23) 



r=2 



with the new coefficients I3\ , e\ and 71 completely determined by the known 



coefficients sj^*'' and Ef' . As the new functions ■d'^x^ (p) and dlj''' (p) are total deriva- 
tives in ^'(p), this equation allows us to calculate ki and Fi. It is obvious that this 
method works for higher g as well. 



3.6. The Iterative Procedure 
With the help of Maple the zero mode coefhcient ki of vi{p) becomes 



lis. 



11 Ao 



Kl 



5A2J2 , 5S2M2 5S2M3 



led^Mi^ 

, 5A2J3 



16d3Ji Sd^Ji-" Sd'^Mi^ 16dMi' 

S2 , A2 ^2J2 



A2M2 



ledJi^ i6d3JiMi md^JiMi ledVi^Afi md^Mi^Ji 
3S2M2' 3A2J2' , 3A3J2 3S3M2 , 5S4 



8 dMi 

5A4 
16dJ^ ' 

3M2^ 

'8¥m7 
11 



8dJi 
5 S3 



idJi^ 8dMi^ 



5 M2 5 J2 



16dA/i' 



^ l8d2Ji 



5M3 



5 J3 



',d^Mi^ 
11 



\d3Ji 



M2 



16d4j/ led^M^^ led^JiMi 



IBd^Mi'^ ledVi^ 



J: 



16d3J 



I'Mi } 



(3.6.24) 



and the doubly fermionic part of Fi is constructed as well. The result for the full 
free energy at genus 1 then reads 



— In Ml ^- In Ji 

12 12 



-(2i-Ai){ 



11 S2 



16^4 Ml 2 
5S2M3 , 5A2J3 



11 A2 

16d4 
^2 J2 



5A2J2 , 5S2M2 



Sd^Ji^ 8d3Mi' 

A2 Af2 



16rf2Mi3 md^Ji"^ Wd^JiMi md^Mi^Ji 



3S2M2' 3A2J2' 



3A3J2 3S3M2 



8^2 Ml"* 

5A4 
16d2 Ji2 

S2 S3 



8d2 Ji^ 

S3 
2d^Mi^ 
A2 A3 



+ 



5S4 



8d2Ji3 8d2Mi3 16(i2Mi2 
A3 , 3S2 3A2 



2d3Ji2 l&d^JiMi l&d^J- 
S2 A2 



— I 

iMi J 



8d2Mi2 8^27^2 Ad^JiMi 



(3.6.25) 



The above results hold true for generic bosonic and fermionic potentials. Note that 
for symmetric bosonic potentials, i.e. x = —y, and generic fermionic potentials one 
has 

Mfc = (-)'+' Jk. (3.6.26) 
easily seen froiH-bhe definitions of eq. (3.5.7). If one takes the fermionic potential 
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to be symmetric as well one finds 



fc > 2, 



(3.6.27) 



from the defining equation (3.5.8). Then one checks = o and this continues 

to hold true for higher genera as well, seen from an analysis of the odd superloop 
equation in the totally symmetric case. 

The Iteration for Ug{p) 

The remaining part of the loop correlator W{ \ p) is now easily derived from Fg 
by applying the loop insertion operator S/dV{p) to its doubly fermionic part. For 
genus 1 the result is 



ui(p) = ^(l«x«(p) + 5«*W(p)), 



(3.6.28) 



where 

r(4) 



A\^> = - 



35 (Si-Ai)S2 



A 



(3) _ 



(Si-Ai 



15 2, 



45 S2 



A 



(2) 



21 (Si -Ai) S3 3 (Si 



2,2 (PMi"^ 
Ai)S2Af2' 



25 S2 Af2 _ 
32d2M,3 32^3 J. 

3(Si-Ai)S2J2 



5A2 -1 
d? Ji Ml i 



"2 "3 



IQcPMi^ 
3 A2 S2 , 



8d2A/i'^ 
3(Si - Ai)S3A/2 



32d3Ji^Mi 8d2A/i^ 
15(Si~Ai)S4 S3S2 



4d^Mi^ IQdPMi^ Ad^Mi^ 

A2(Si-Ai)A/2 A2(Si-Ai) 



1« = 



8 JiArid3 

51(Si-Ai)S2 

32d^Mi^ 32 Jid^Mi^ AJid^Mi 

33 (Si-Ai)S2Af2 ^ 5(Si-Ai)S2A/3 _ 9 (Si 
32d3Mi^ 

33(Si-Ai)S2Af2 (Si 
32d4Mi^ 

3(Si-Ai)S2J2^ 

5(Si~Ai)S2 7 (Si 
16d5JiA/i 32d4Ji' 



Ai)S2 



lad^Mi^ 32Ji'a*Afi 

■Ai)S2Af2 35 (Si-Ai)S5 



32d4j^Mi^ 



32^2^1^ 

15(Si-Ai)S4A/2 , 45(Si-Ai)S4 



16 d^Mi^ 
Ai)A3 



So A3 



32d3Mi^ 

^3 "2 



Wd^Ji^ 16d^Mi' 
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51(Si-Ai)S2 _ 51(Si~Ai)S3 _ S2 S3 , 5 S2 S4 
32d^Mi^ 32d4Mi2 



O-L y^l — ^2 ^J- V^l ~ ^'-i; ^3 ^2^3 o^2^i 

32 d^Mi^ 32^4 Ml 2 8d3Mi^ IGd'^Mi^ 

A2 S2 _ 7A2 (Si- Ai) J2 3A2 S3 7A2 (Si - I 

"2JiAfid4 32Ji^d^ SdUiMi IGJi^d^^ 
, 5(Si-Ai)S: 



16 Ji^d^ 



7 A2 (Si -11)1 

SdUiMi 16Ji2d5l=l 

■.2 J2 ll(Si - Ai)S2 _ 5(Si-Ai)S2 Af3 

,0,r- -'3 



9(Si-Ai)S3 A/2' 



32 Ji^ds 



16 d3Mi 



3(Si-Ai)S3J2 9(Si-Ai)S3 

32d3Ji2Mi 32JiAfic«4 

9(Si-Ai)S3Af2 ^ 5(Si-Ai)S2J3 

8d3Mi3 32d3j^3 



16d2Mi4 
^ 15(Si-Ai)S3Af3 
32 d2Mi^ 

S2 / 11 A2 5 A2J2 5 A2J3 A2 
~2d \lQd^Jx^ ^ 8dUi^ ^ 16dJi^ ^ 16 dMiMi 

A2M2 _ 3A2J2^ ^ 3A3J2 _ 3 S3A/2 ^ 5S4 
16^2 Ml Vi 8dJi'' 8dJi^ 8dA/i^ 16dMi^ 

5A4 5S3 5A3 \ ^ 3(Si- Ai)S2M2" 

16dJi^ 8d2A/i2 8d2JiV 8^3 Mi'' 

(Si-Ai)S2J2 ^ A2(Si-Ai)Af2 ^ llA2(Si-Ai) 
Sd-^Ji^Mi led'^JiMi^ 32Jid5Mi 

f" - A 1 1 3 J2^ _ 3M2^ _ 5Af2 _ 5 J2 

^^^^ ''Xsd^ji^ sdHii^ sdnii^ 8d3j,3 
5 A/3 5J3 11 11 

-I i 2 1 

16d2A//i3 16d2j^3 16d4j^2 IQd^Mi^ 

16d^JiMi^ IQd^Ji^Mi / 



and the analogous expressions for the o'f^ obtained 
M ^ J, E ^ A and d -d. 



(3.6.29) 

from the above by replac: 



;ing 



General Structure of Ug, Vg, Ug and Fg 



In the following subsection we deduce the number of moments and basis functions 
the free energy and the superloop correlators at genus g depend on. 

Ambj0rn et al. [8, 39] have shown that the free energy of the hermitian matrix 
model depends on 2(3(7 ~ 2) moments. This directly translates to F^"". Similarly as 
Ug = 5Fg°^/5V{p) and with eq. (3.6.11) we see that Ug contains 2{3g — 1) bosonic 
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moments and basis functions up to order (Sg — 1), i.e. 



3g-l 



(3.6.30) 



fc=i 



For the structure of Vg consider the leading-order poles on the right hand side of 
eq. (3.6.2). Label this order by Ug, then with eqs. (3.6.10), (3.6.11) and (3.6.13) the 
three terms on the right hand side of eq. (3.6.2) give rise to the following poles of 
leading order 

W o ugip):{3g-l) + l 

Vg' Ug-g' ! + 3 (ff ^ ) - 1 ^ + 1 



From the above we deduce that Ug — 3g, and therefore 



(3.6.31) 



3s 



(3.6.32) 



k=l 



As the highest bosonic moments in Vg come from the highest-order basis functions, 
we see that Vg depends on 2(35) bosonic moments. To find the dcpeadcncc on 



Jiavc 



the number of fcrmionic moments recall eq. (3.6.7). In order for Vg 
leading contribution of x'^^^Hp) the fermionic part of the free energy P'f"'^ must 



contain E^g+i. We are thus led to the conclusion that Fg'^""^ and Vg both depend on 
2(35 + 1) fermionic moments ^. As the application of the loop insertion operator 
6/5'^{p) does not change the number of bosonic moments, Fg'^"'^ must contain 2(3 g) 
bosonic moments. 

Knowing the structure of F^^""^ then tells us with eqs. (3.6.10) and (3.6.11) that 
Ug depends on 2(3 g + 1) bosonic and 2(3 5 + 2) fcrmionic moments. For genus g it 
reads 



35+1 



(3.6.33) 



k=l 



^In some sense this is counterintuitive: In ref. [8] it was shown that there is a relation of the 
dependence on the number of moments to the dimension of the moduH space of a Riemann surface 
of genus g, which is 3g — 3. So one could have speculated that something similar hol ds true for 
the fermionic moments and the number of fermionic moduli of a supcr-Riemann st rfa ;c, which is 
2g — 2. However, the discovered slope of 3 in g seems to contradict such an intorpi eta ;ion. 
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The genus g contribution to the ( |s)-superloop correlator Wg{ \p, . . . ,p) will then 
depend on 2(3 g + s) bosonic and 2(3 5 + s + 1) fermionic moments. Similarly the 
genus g (l|s)-superloop correlator Wg{p\p, ■ ■ ■ ,p) is a function of 2(3.g + s) bosonic 
and 2(3 g + s + 1) fermionic moments. |^ 

This concludes our analysis of the iterative process away from the doubleCEEB-ling 
limit. 



7. The Double Scaling Limit 

In view of our lack of understanding the supereigenvalue model on a geometrical 
basis as a discretization of super-Riemann surfaces it appears quite ambitious to 
speak of a continuum limit. However, due to the striking similarities to the Hermi- 
tian matrix model with its well understood continuum limit (discussed in chapter I 
section 2) it is natural to export these techniques to the supersymmetric case. We 
shall see that again one finds critical values for the 1 
allowing one to perform a double scaling limit. 



CO ipling constants gk and ^fc+1/2 



Similar to the situation in the Hermitian matrix model the "naive" — > 00 limit 
of the supereigenvalue model is unsatisfactory as it leaves us only with the planar 



contributions, easily seen from eqs. (3.4.8) and (3 



). Crucial for the double scal- 



ing limit is the observation that there exists a sr bs )ace in the space of couplings 
{gk, ^fc+1/2} where all higher genus contributions ttUjie free energy Fg diverge. This 
enables us to take the double scaling limit, where one simultaneously approaches 
the critical subspace of couplings fud takes N ^ 00, giving contributions to the free 



tl e Hermitian matrix model Kazakov [12] multi- 
a zeros of the eigenvalue densities accumulating 



energy from all genera. Just as ii 
critical points appear, related to extr 

at one endpoint of their support. 

For the formulation of this limit the moment technique developed in the above 
turns out to be extremely useful. One can then determine which terms in the explicit 
solutions for Fg, W{p\ \ ) and W{ \ p) calculated in the previous section contribute 
to the double scaling limit. However, the calculation of these quantities away from 
the double scaling limit is rather time cogsgiming, especially as a large number of 
terms turn out to be irrelevant in this limit. In this section wc will therefore develop 
a procedure which directly produces only the terms relevant in the double scaling 
limit. 

The first thing to do is to study the scaling of the moments and basis functions 
by approaching the critical point. 

Multicritical Points 



The analysis of the scaling behaviour for the bosonic quantities was carried out by 
Ambj0rn et al. [8, 40] in the framework of the hermitian matrix model. Consider 
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the case of generic, i.e. non-symmetric, potentials V{p) and ^'(p). The m'th multi- 
critical point is reached when the eigenvalue density {uo{p) — V'{p)) of eq. (3.5.22), 
which under normal circumstances vanishes as a square root on both ends of its 
support, acquires (m — 1) additional zeros at one end of the cut, say x. Alterna- 
tively one would define the critical point by a diverging free energy. A glance at 
our genus one result of eq. (3.6.25) tells us that Fi diverges for vanishing Mi (or 
Ji). The only moments appearing in the denominators of higher genus Fg^s will be 
Ml and Ji as well. Hence the condition for being at an m'th multicritical point is 
taken to be 

M^ = M^ = ---M^^_,=0, Aq^O, k>m, Jf ^ 0, / > 1, (3.7.1) 

defining a critical subspace in the space of bosonic couplings gk- Denote by g'j: a 

distribution is confined 
int the conditions of eq. 



particular point in this subspace for which the eigenvalue 



to the interval [xc, yd- If we now move away from this 
(3.7.1) will no longer be fulfilled and the cut will move to the interval [x, y]. Assume 
we control this movement by the parameter a and set [8] 

a; = a;e-aAi/" 

p = Xc + aTT (3.7.2) 

The scaling of p and the introduction of its scaling variable tt is necessary in order 
to speak of the double scaling limit of the superloop correlators. A plays the role of 
the cosmological constant. We will now deduce the scaling of y by further assuming 
that the critical subspace {g'j.} is reached as 

gk^g- gl (3.7.3) 

where g is a function of a to be determined. Imposing the boundary conditions 
(3.5.3) yields 



y — yc ^ a™, and g — I ^ a" 



(3.7.4) 



This may be seen as follows. Consider the first boundary condition (3.5.3) using J/ = j/c + ^U, 
eqs. (3.7.2) and (3.7.3) and the fact that the for i S [l,m — 1] vanish lead to 

0- '^'^ S^c'H 



^2n (u;- -|- a Ai/'")i/2 {u) - yc - Aj/)i/2 

oo , , . oo 



^ i— i— 
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oo 

= g ^ c, Mf (aAi/™)' +gc^ Ay 

i — m 

+0{a-Ay), 

from which one deduces Ay ~ a™ . The computation of g — 1 ~ a™ goes along the same lines 
starting from the second boundary condition (3.5.3). 

Knowing this one easily computes the m'th multicritical scaling behaviour of the 
bosonic moments 



Mfe-a™-^ /cG[l,m-l], 
while the higher Af -moments and the J-moments do not scale. 



(3.7.5) 



See this by direct computation: 

Ml, 



p 2n (<.^-Xc + aAl/'")'=+l/2 (cj-yc-a™A)l/2 



v;'(<^) 



(aAV'")' ^^V^ (0) (a^Ky 



{UJ - XcY 



j =0 



where for £ [1, m — 1] the leading nonvanishing term is a™ ^ (from i = m — k and j = in the 
sums). I I 

Moreover the functions (^i"'' (p) and (j}^^ (p) are found to behave like 

4")(p)^a-"-l/^ 4")(p)^a-i/2, (3.7.6) 
from which one proves the scaling behaviour of the basis functions 

x("Hp) ~ *(")(p)-a-i/2^ (3.7.7) 

by induction. 

Let us now turn to the scaling of the fermionic moments and A^. Similar to 
the bosonic case the function {vq{p) — ^'(p)) of eq. (3.5.24) usually vanishes at the 
endpoints of the cut like a square root. We will fine tune the coupling constants 
Cfc+i/2 in such a manner that (n — 1) extra zeros accumulata-at x, i.e. 



"2 



0, 



S^T^O, k>n, Ai^O,^ 



> 2, 



(3.7.8) 



where 5J: = .^.^[ccc, 2/c, ^^+1/2]- addition the analysis of tlpe-eolution away from 
the scaliiiff limit tells us that the moments St— and Ai will alvays appear in the 
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3. 7. The Double Scaling Limit 



combination (^i — Ai) ^. This suggests to impose the constraint — = on 



these moments. As there is no analogue to the bcurdary conditions (3.5.3) for 
the fermionic quantities we are free to choose the scaliag of the couphng constants 
^fc+i/2- We set 



+ 1/2 



,1/2 



fc+1/2' 



derives (Si — Ai) 



(3.7.9) 

„n-l/2 



and will comment on this choice later on. From this c ne 

and Sfe ~ gn-fc+i/2 ^ g [2, n — 1]. All other fermionic moments scale uniformly 
with a^/^. So far the fermionic scaling is completely independent of the bosonic 
scaling, governed by the integer n. We shall, however, introduce the requirement 
that the scaling part of the bosonic one-superloop correlator Wo{ \ p) of eq. (3.5.21) 
scales uniformly, i.e. we require {uo{p) — V'{p)) and uo{p) to scale in the same way 
[31]. As {uo{p) — V'{p)) ^ a™^^/^ we arrive at the following condition on n 



1. 



(3.7.lmj 



And therefore 



Si-Ai 



1+1/2 



k> 



m-k+3/2 



(3.7.11) 



and a — s- but keeping 



The double scaling limit is now defined by letting N - 
the string coupling constant a — a^^™"^ iV~^ fixed. ^ 

The Iteration for Ug{p) and Vg{p) 

By making use of the above scaling properties of the moments and basis functions 
we may now develop the iterative procedure which allows us to calculate directly 
the double scaling relevant versions of Ug (p) and Vg (p) . Recall the iterative scheme 
of the previous sections: By eq. (3.5.13) every Ug and Vg may be written as a linear 
combination of basis functions x^"-* and \I>("), where the coefficients of this expansion 
are simply read off the poles at x and y of the right hand sides of the superloop 

'^The only exception is the planar zero mode x of cq. (3.5.2Cl^^owever, as seen from eq. (3.5.24) 
vo{p) effectively does not depend on this combination. 

*Let us now comment on the choice of eq. (3.7.9). At first sight one might have expected 
a scaling like Jfc_|_i/2 = ~^ °{'^)] ^^+1/2' '^^^^ turns out to be inconsistent because then the 
condition (3.7.10) demands n to be half— integer which it can not be. If one takes the more 
general ansatz = ?fc+i/2 condition of uniform scaling of W( | p) yields the allowed 

sequence {n, 1} = {m + 1, 1/2}, {m, 3/2}, {m — 1, 5/2}, . . .. The scaling of the lowest moments in 
eq. (3.7.11) remains unchanged, the uniform scaling however already starts with H,i scaling like a' 
and thus simply reduces the number of double scaling relevant terms. 



equations (3.6.1) and (3.6.2) after a partial fraction decomposition. To optimize the 
procedure to only produce terms which are relevant in the double scaling limit we 
have to analyze the operators appearing on the right hand sides of the superloop 
equations, i.e. the superloop insertion operators S/5V{p) and 6/6'i{p) as well as the 
operator W. 

From the point of view of the a ^ limit the effect of a given operator in S/SV{p) 
acting on an expression which scales with a to some power is to lower this power 
by a certain a mou nt. Carefully examining each term in 5/dV{p) shows that a is 
maximally lowe'rsi by a power of (m + 3/2). All operators which do not lower a 
by this amount are subdominant in the scaling limit and may be neglected. The 
outcome of this analysis for S/6V{p) is 



5Vip)^ 



E 



SMk 



^ - 5Vip)^ SMk 
J(Si-Ai) 



Sx d 
6V{p) 

6 



OO r-,— , 



SV{p), (5(Si-Ai) 



E 

k=2 



6V{p)x 5^k 



where 



SMk 
SV{p), 



+ (fc+-) 



M, 



k+l 



Ml 



(3.7.12) 



(3.7.13) 



S^ 

SV{p) 



Ml 



D 



SiEi - Ai) _ 1 S2 



SVip)^ 
SVip)., 



2 Ml 



2 



Mr 



and 



(3.7.14) 
(3.7.15) 



4")(P) = (P-X) — V2rf^-l/2^ 

with dc = Xc — Vc- Repeating this analysis for the fermionic superloop insertion 
operator S/S'^iji) shows that here a is maximally lowered by a power of (m + 1) 
and the relevant contributions are 



^(Si-Ai) 



S'i!{p)x 



S^{p).c S{Si - Ai) ^ S^{p)x SEk 



E 



SE. 



(3.7.16) 



'Here the subscript x indicates that the critical behaviour is associated with the eildpDint x 



jidp 



Ill: The Supereigenvalue Model 



with 



(3.7.17) 



Finally we state the double scaling version of the operator ^ snOng on the function 



o 



n „ 



1 



fc=i 



(Si-Ai) 

2dr 



1 



(p-x) 



n+1 ' 



(3.7.18) 



Here the operator is seen to increase the power of a of the expression it acts on 
by a power of m. 

We are now in a position to calculate the double scaling limit of the right hand 
sides of the loop equations (3.6.1) and (3.6.2) for Ug{p) and Vg{p) provided we 
know the double scaled versions of Wi(j3), . . . ,Ug_i(j)) and Vi(p), . . . ,Wg_i(p). As 
all the y dependence has disappeared we do not have to perform a partial fraction 
decomposition of the result. Moreover no Jk and Afc dependent terms will contribute 
if we do not start out with any and we do not. 

The starting point of the iterative procedure are of course the genus correlators. 
Keeping only the double scahng relevant parts of eqs. (3.6.5) and (3.6.8) one has 



d dvoiq) 



1 



dq S'^ip)^ 
Similarly eq. (3.6.12) turns into 



\p=q 



Suq{p) _ 1 

6V{p)^ " 8 (p - x)2 



Svojp) 
SV{p)^ 



(Si-Ai) 



4 dc Ml 



{P - 



^2 



4 dc Ml 



{P - 



(3.7.19) 



(3.7.20) 



The higher genus correlators Ug (p) and Vg (p) are expressed as linear combinations 
of the basis function x'^^^Kp) take the general form 

3s- 1 
k=l 



3g 



,{p)=Y,bI,'^X^'Hp) + ng4,("\p), 



(3.7.21) 
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where Hg is the zero mode coefficient not determined by the first two superloop 
equations (3.6.1) and (3.6.2). Note that the Ag'^' coefficients should up to a factor 
of two be identical to those of the double scaled Hermitian matrix model obtained 
in ref. [8]. We have calculated the Ag''^ and Bg'''' coefficients in the double-scaling 
limit for 5 = 1, 2 and 3 with the aid of Maple. The results of Ambj0rn et al. [8] for 

the Ag''^ coefficients could be reproduced. 

Before we state our explicit results let us turn to the general scaling behaviour of 
the one-superloop correlators 



W^<,(|p)^a(i-2<,)(™+i/2)-i^ l^<,(p|)^a(i-2<,)(™+i/2)-i/2^ (3 7 22) 



which one proves by induction. As shown in ref. [8] the coefficients Ag'''^ of eq. 
(3.7.21) take the form 



Ma, ■ . . Mg^ 

Ml" 4^-1 ' 



(3.7.23) 



where the brackets stand fot^tional numbers and where aj,a and s are subject 
to the constraints 



a = 2g 



,nd 



i=i 



a, 



1) = 3.g - fc - 1 



(3.7.24) 



with aj e [2, 3g — 1]. Explicitly up to genus two one finds 
1) 



= 



4(2) _ i 



A^i^ = 



A 



(3) 



49 Mo^ 



49 M2 
64 dcMi= 



105 



128 da Ml 



2 ■ 



5M. 



128 dc Ml'' 16 dc Ml' 



(3.7.25) 



The terms above are only potentially relevant, as for the m'th multicritical point 
all terms containing Mk with k > m arerSHibleading in the double scaling limit. 

(k) 

Similarly the structure of the coefficieitf Bg may be determined from the iter- 
ative procedure and is seen to be ' — 



a,, 13 



■,01S,P \ Oi)g^k 



. M„^ Ef3 



Ml" dc 



(3.7.26) 
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where the brackets denote rational numbers and where we write Si for (Si — Ai). 
One shows that the a, aj, (3 and s obey the conditions 

s 

a = 25 + 5-1, and ^(a^ - 1) = 3.g + 1 - /3 - fc (3.7.27) 

with aj e [2, 3.g] and /? e [1, 3(7]. For the zero mode coefficient Kg the general struc- 
ture is given by the same expansion as eq. (3.7.26) with k = Q. The conditions on 
a, Uj, (3 and s then read 

s 

a = 2.g + s, and ^(a^ - 1) = 3.g + 1 - /3 (3.7.28) 

where aj € [2, 3(7] and (3 G [1, 3^ + 1]. 

The explicit results for the Bg'"'' coefficients ioi g — 1 and g — 2 are given by 

^(1) _ ^3 , M2E2 



1 SAhdc SAfi^d,' 

(2) _ M2(Si - Ai) S2 (3) _ 5(Si-Ai) 

' leMi^de 8Mic;,' 1 16 Mid, ■ 

^(1) ^ 203M2S5 _ USMs'Sa _ 105 Sg _^ 63M2'S3 
^ 128 4^ 128Mi^(fc^ 128 4^ Mi^ 32Mi'^dc^ 



B2 



(2) 



105M4S3 


145 M3S4 105 M5S2 


77M4M2S2 


128 4^Mi'* 


128 4^Afi'* 128dc^Mi^ 


32Mi^4^ 


87M3M2S3 


75M3M2^S2 63M2^S2 


63M2^S4 


32 Mi^dc^ 


WMi^dc^ 32Mi^(fc^ 


32 Ml 5 4^' 


21 M2^S2 


2IM2S4 77M3M2S2 


35 M4S2 


GAMi^dc^ 


64 4^A^i* 128Mi^(fc^ 


USdc^Mi'^ 


75M3M2^(S 


i-Ai) 105M5(Si-Ai) 


63Af2*(Si - 



32Mi^dc^ 256 4^Mi^ 64Mi^4^ 

35 S5 145Af3'(Si- Ai) , 21M2'S3 35Af3S3 



B. 



USMi^dc^ 256Mi^dc^ QAMi^dc^ 128 d^^Mi^ 

_77MiM2{Ei - Ai) 
64M?d? 

(3)_ 599M3M2(Si - Ai) 105M23(Si- Ai) 5M3S2 



^ 128Mi^dc^ 32Mi^dc^ Wdc^Mi* 
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385M4(Si - Ai) _^ 21M2^S2 , 7Af2S3 



B2 



B2 



256 4^7^1"* GAMi^dc^ 128 4^Mi*' 

(4) 357 1/2' (Si -Ai) ^ 875M3(Si - Ai) _ 35 S3 _ 63M2S2 
64Mi^d,^ 256 dJ 

(5) _ 105 S2 1617Af2(Si- A] 



Wi^ 128 Mi-'dc 128 d/Mi^ 



128Mi^dc 256(ic Ml 



^(e)^_1155(^i_Ai)_ (3^^ 
^ 256Mi^(fc 

Note that the terms listed above are only potentially relevant, depending on which 
multi-critical model one wishes to consider. For an m'th multi-critical model all 
terms containing Mfc, fc > m, or S;, I > m + 1, vanish in the double-scaling limit. 
We remind the reader that we assumed to have a non-symmetric potential and that 
the critical behaviour was associated with the endpoint x. In the case where the 
critical behaviour is associated with the endpoint y all formulas in this section still 
hold provided dc is replaced by —dc, Mf; by J^, S^. by A^, and x by y. 

The Iteration for Fg, Kg and Ug{p) 

Having computed Ug{p) and Vg{p) (up to the zero mode) we may now proceed 
to calculate the free energy Fg and the zero mode coefficient Kg. This is done 
by rewriting Ug and Vg as total derivatives in the superloop insertion operators 
5/SV{p) and S/S'if{p) respectively, yielding the bosonic and doubly fermionic parts 
of Fg as well as Kg. The procedure to compute the bosonic part of the free energy 
Fg°^ works just as in the Hermitian matrix model described in ref. [8]. We have 

From eq. (3.7.22) we see that Fg scales as 

Fg = Fl^°'' + Fgf™ _ a(2-23)(m+i/2)^ (3.7.30) 

just as the hermitian matrix model at its m'th multicritical point. 

To find the bosonic piece of the free energy Fg one uses the following recursive 
rewriting of x*^"Hp) ^ total derivative 



n-1 



Ml V 2n - 1 6V{p) 



k=2 



One does not need any expression for x^^\p) ^ ^s^"* = in eq. (3.7.21) for all g >2 
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To obtain the doubly fermionic part F^"^ simply rewrite the basis functions 

of Vg{p) appearing in (3.7.21) in terms of the functions 4>''J^\p) which by eq. 
(3.7.17) are nothing but total derivatives in 6/6'i/{p). Doing this for Vg{p) allows 
one to directly deduce the form of and Kg. ^° 

The explicit results for genus one are 

F^^ = -i-lnMi, 



fcrin /<— ' 



(Si - Ai 



)(- 



5S4 



3M2S3 



5M3S2 3M2^S2' 



and 



Ki = 



—2 ^3 

5S4 



3 M2 S3 5 M3 S2 



+ 



3M22S2 



(3.7.32) 



(3.7.33) 



WMi^dc SMi'dc 16 Mi^dc SMi^dc 
These expressions may of course be alternatively obtained by taking the double 
scaling limit of the genus one results computed away from the scaling limit, except 
for the scaling violating term Fi°^ at genus one. This term in eq. (3.7.32) was 
selected by producing the double scaling version of ui{p). 

Before presenting explicit results for genus two, let us describe the general struc- 
ture of the free energy at genus g. For the bosonic piece we have [8] 



E 



ai,. . .,as 



a, 



. . . Af„^ 



5>2, 



M?'dc^-^ 

Qj > 1 1 1- 

where the brackets denote rational numbers. The a-,-, a and s obey 



(3.7.34) 



2g - 2 + s 



and 



E("J 



1' 



3,9 - 3, 



with Uj e [2,351-2]. 

Similarly the general form of F^'^'' 



rfcrm 



, as,/3i,/32 I a) 



for (? > 1 reads 

^/3i 



.Mr, 



Ml°'d, 



9+1 



(3.7.35) 



(3.7.36) 



^"Actually the presentation here is somewhat misleading. In order to compute the Ug{p) and 
Vg(p) by iteration in the way outlined in the previous subsection one needs to know the full 
Vg—i{p) including the zero mode coefficients ki, . . . , Kq-i. In practice one hence com- 



putes all quantities at genus g, i.e. Ug, Vg, Kg and Fg, before proceed 



ng 



to genus (g + 1). 
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Dou 



)7e 



Scaling Limit 



where we write .^1 for (^1 — Ai) and where the brackets eenote rational numbers 
as before. The aj, Pi, a and s are subject to the constraints 



a 



2g + s, and 



where aj G [2, 3g] and Pi £ [1, 3^ 
The results for genus two read 



29M2M3 21 M2^ 

^2 = 



E 



(aj-l) = 3.g + 2-/3i-/32, 



(3.7.37) 
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4 Ml 
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693 Ee M2 
64MA^ 
_ II55S3M2M4 
64Afi^4^ 



315 S3 Ms 2175S3Af3^ 675S3M2^M3 1785S2M3M4 



64Mi^dc 



256 Ml 



16Mi^4 



128 Mi^dc 



495S2M2^M3 2205S2M22M4 1323 S3 M2'' 



8175S2M3^M2 



8Mi*dc^ 
1155 S2M6 
256M7d7 
105 S2 S4M2 
64Mi^4^ 



64Mi^dc^ 64Mi^4^ 256Mi^4^ 

2IS2SSM2 , IG5S2S6 I45S2S4M3 



+ 



16 Ml 

21 S3 S4 M2 

64Mi^4^ 128Mi''4^ 32Mi^dc 



128 MfdJ^ 
35 S3 Ss 



128Mi^4^ 
63 S2 S3M2^ 



35 S2 S3 M4 _^ 195 S2 S3 M2 M3 
32Mi^4^ 



64 Ml'' 4 



(3.7.38) 



Of course so far we have determined only the coefhcients Fg of the genus expansion 
of the free energy (cf. eq. (3.4.9)). For an m'th multicritical model the relevant 
expansion parameter in the double scaling limit is the string coupling constant 
a — a"^™"^ N~'^. If we introduce the bosonic and fermionic scaling moments /ifc 
and Tk by (cf. eqs. (3.7.5) and (3.7.11)) 

Mk^a"'-'' fik, fce[l,m]. 



(Si - Ai) = a 



= „™+i/2 



Tl, 



S,=a"-'+3/2r,, 



I e [2,771 + 1], 



(3.7.39) 
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we get by replacing Mk by /i^ and setting Mk equal to zero for A; > m as well as 
replacing S; by ti and setting S; to zero for Z > m + 1 in the formulas above exactly 
the coefficients of the expansion in the string coupling constant. Needless to say 
that these results apply for non-symmetric potentials as well, where the critical 
behaviour is associated with the endpoint y by performing the usual replacements. 
For the zero mode coefficient at genus two we find 



K2 



495 SaMaMz^ BSSaAfsAfa 315S3A//5 1323 SgMa^ 
GTSSgAfaMs^ 2175S3Af32 8175 SaAfa^Afa 385S4Af4 



ledc^Afi^ 256 4^Afi^ 256 dc^Mi'^ 64 dc^Mi^ 

2205S2A//4Af2^ 375S4Af3Af2 1155S2Af6 1785S2Af4Ar3 
64 4^Afi^ 16dc^Mi^ 256d7M^ 128 d^^ Mi 

693S6Af2 315S4A/2^ 525S5Af2^ IOI5S5A//3 



6 



64 4^Afi^ 16 4^A/i^ 32dc^Mi^ 128dc^Mi^ 
1323S2Af2' ^ 1155S3Af4Af2 ^ 1155^7 
64 4^A/i^ 64 4^A//i^ 256dc^Afi^' 

We obtained these results with the aid of a Maple program which performs the 
iteration up to arbitrary genus. In practice the expressions become quite lengthy, 
e.g. consists of 114 terms. 

For the knowledge of the full Wg{ \ p) it remains to compute Ug{p). This is of 
course done by applying S/dV{p) to . The general structure of Ug{p) turns 

out to be 

3g+l 

(3.7.41) 

k=l 



where 



AW= K,...,a.,A,/32|a),,. ^"-^^f;;:f (3.7.42) 



underlying the conditions 

s 

a = s + 3.g, and ^(oj - 1) = 4 + 3.g - A: - /3i - /32, (3.7.43) 

with aj e [2, 3g + 1] and (3i £ [1, 3^ + 2]. Due to space let us only state the genus 
one results 

-(1) _ 5S2 S4 ^ 3Af2S2S3 _ 35 (Si - Ai)S5 
^ ~32 4^Afi2 16 4^Afi^ 32 4^Afi^ 
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15Af3(Si 



Ai)S3 



15A/2(5i - Ai)S4 9Af2'(Si - Ai)S3 



A 



(2) 
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.n 



16 4 Afi^ 



^(3)— _15(^l-Ai)^3 
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2/^ 



\dc^Mi 



Adc^Mi^ 
Ai)S2 



A 



(4) 



16 4^Afi^ 



35 (Si -Ai) 52 
32 dc^Mi^ 



Sdc^Mi^ 
25Af2(Si - Ai)S2 
32 4^A/i^ 



(3.7.44) 



This concludes our analysis of the double-scaling limit for generic potentials. 

The case of symmetric bosonic and generic fermionic potentials was considered 
in ref. [30] where a doubling of degrees of freedom was observed for genus zero. 
We then have the independent set of moments {Mk,^k, ^k}- With the methods 
presented in this paper one can see that this holds for higher genera as well, i.e. the 
free energy here takes the form 



^symm _ 2 ^bo 



F; 



ifcrm,: 



pfcrmiy- 
9 



(3.7.45) 



where Fg°^ denotes the bosonic part of the free energy for generic potentials, 
^form,x g^^^ ^fcim,y=-x (jgnQ^g ^hc doubly fcmiionic parts in the generic case where 
the critical behaviour is associated with the endpoints x and y = —x respectively 
(cf. eq. (3.7.36)). If one chooses to take the fermionic potential to be symmetric 
as well the doubly fermionic part of the free energy will vanish, as we have already 
seen on the discrete level in section 6. 



8. Identification of the Model 

In order to identify the continuum theory described by the supereigenvalue model we 
will have to discuss some results obtained in super-Liouville theory. This subsection 
will by no means give a review on this subject, we simply highlight some formulas 
relevant for the interpretation of our results. For more details see refs. [41, 42, 4]. 

In the following we consider N = 1 superconformal field theories with c = d cou- 
pled to 2d supergravity. For a base manifold of fixed genus g the partition function 
is given by 

^9=E JWEtj][VX']e-''i''^^\ (3.8.1) 
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where the sum runs over the spin structures. Here E^^ denotes the super-zweibein, 
the free, conformal superfield and the matter action reads 

S'[X,S] = i / d^zd^OD^X^ D°'X^, (3.8.2) 

using the superspace notation of ref. [41]. Here Da are superdifferentials and 
/ = 1, . . . ,d. The action and the measures DE^j and VX' are invariant under 
superdiffeomorphisms and local Lorentz transformations. These have to be divided 
out of the path integral (3.8.1) which is symbolized by square brackets around the 
volume elements. The action (3.8.2) is also invariant under the group of local super- 
conformal transformations while the measure is not. Thus a superconformal gauge 
fixing to the reference frame E — e"''^ E will yield a Jacobian which is assumed to 
be [41] 

J(0,^) =exp(-5sL[0,^]). (3.8.3) 
The action functional in the exponent is the super-Liouville action 

Ssl[^,E]^^ J d^zd^9E{^Da^D^(P~^ ^ Y(f> + K a.in e'' * ) (3.8.4) 

where (p is the Liouville superfield and Y the curvature superfield. The cosmological 
constant k couples to the operator of lowest-dimension in the Neveu-Schwarz sector 
of the matter theory denoted by Omin [43]. For unitary theories Omin is unity. 
Moreover a is a c dependent constant. The precise form of the matter action (3.8.2) 
is not of importance, all that is needed is that one has some superconformal field 
theory with central charge c. One can now study the partition function (3.8.1) in 
the form 

Zg= [ dLZg{L), ^ (3.8.5) 

Jo 



where L is a characteristic length scale of the super-Riemann surface 
partition function for the super-random surface of fixed length is defined to be 

Zg{L) = Y, J T)X'D<f)e-^si.-SM5(^j Szd^OE-L 



(3J8.6) 



where the symbol DX denotes the integration over all the super-matter, super-ghost 
and super-moduli contributions, contained in the matter action Sm- Similarly to 

^'^Note that L has the dimension of a length. In contrast to the bosonic random surface which 
has an intrinsic area defined by J d?z the analogous object for a super— random surface is given 
by J d?zd?6 E = L. Because of the Grassmann integrations the dimension of L is lowered by 1 
compared to the bosonic case. 
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3.8. IdentiGcation of the Model 
the bosonic case the scaling behaviour of Zg{L) may be computed [41] 

Zg{L) = e(«--«) ^ (7str.-2)-l J^^ (3 8 7) 

where and fg are undetermined and L independent. Here the string susceptibility 
7str. of a unitary superconformal field theory takes the form 



7str. 



v/(9-c)(l-B)). 



(3.8.8) 



The minimal iV — 1 superconformal field theories with c < 1 are classified by a pair 
of integers (p, q) and have central charge 



c= 1 - 



)Cw 



2(p-g)^ 

pq 



(3.8.9) 



The unitary theories (with c > 0)L-a¥e given EjZlthe pair {p, q) = {m + 2, m) with 
m > 2. For a non-unitary theory formula (3.8.8) is no longer valid and instead one 
has [43, 4] 

7s.. = 2-li^. (3.8.10) 

p + q-2 

In any case the integral over L in eq. (3.8.5) with the Zg{L) of eq. (3.8.7) may be 
performed to give 



{k - K,)(2-7str.) (1-9) Y[{g _ 1)(2 - ^,,,,)] fg 



(3.8.11) 



note the analogy to the bosonic result of eq. (1.2.11). Guided by the approach in the 
Hermitian matrix model, one now identifies the free energy of the supereigenvalue 
model with the partition feaction of 2d supergravity coupled to superconformal 
matter, i.e. 

N^Fg[gk,^k+i,2]=Zg, (3.8.12) 

and sets {gk — gj:) ^ [n — Kc). Using our result of eq. (3.7.30) and the scaling be- 
haviour of the bosonic couplings in eq. (3.7.4) one finds 



(3.8.13) 



hence we see that 7str. = — 1/m. The computation of a string susceptibility ex- 
ponent is of course not enough to identify the continuum model described by the 
double scaled supereigenvalue model. One has to establish a connection between Li- 
ouville operators to supereigenvalue correlators. The Neveu-Schwarz and Ramond 
sectors of the super-Liouville theory are refiected in the even and odd correlators 
cjftjie supereigenvalue model A^'^) and ■^i'^) respectively. Adequate linear 
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combinations of these correlators in the scahng limit can be identified with super- 
Liouville amplitudes, as was shown by Abdalla and Zadra in ref. [36] by using the 
planar supereigenvalue results of ref. [30]. The outcome is that the double scaled 
supereigenvalue model at its m'th critical point describes the coupling of a minimal 
superconformal field theory of type (p, q) — (4m, 2) to 2d supergravity having the 
central charge 

(2m -1)2 7 22 , , 

d=l-^ L = o,--,--,... m>l. 3.8.14 

m 2 3 

The only unitary model in this series is pure supergravity with m = 1. Note, how- 
ever, that this case must be treated separately, as in our approach m — 1 corre- 
sponds to no restrictions on the moments (cf. eq. (3.7.1)), and thus leads to no 
critical behaviour. This is completely analogous to the bosonic case, where a string 
susceptibility of 7str. = — 1 corresponding to c = —2 may not be reached by the Her- 
mitian one matrix model, but rather through the model studied in chapter II. For 
the higher non-unitary superconformal theories one easily convinces oneself that 
(p, q) = (4m, 2) yields 7str. = — 1/m by using eq. (3.8.10)^ 

9. Conclusions 

We have constructed the supereigenvalue model by imposing the super- Virasoro 
constraints on its partition function. By integrating out its fermionic degrees of 
freedom we saw that the bosonic part of the free energy of the supereigenvalue 
model is simply twice the free energy of the corresponding hermitian matrix model. 
First the model was solved away from the double scaling limit. The super loop 
correlators obey a set of integral equations, the superloop equations. These two 
equations could be split up into a set of four equations, jSorted by their order in 
fermionic coupling constants. By a change of variables froiH-' coupling constants to 
moments we were able to present the planar solution of the superloop equations 
for general potentials in a very compact form. The remarkable structure of the 
superloop equations enabled us to develop an iterative procedure for the calculation 
of higher-genera contributions to the free energy and the superloop correlators. 
Here it proved suSicient to solve the two lowest-order equations at genus g for the 
purely bosonic Ug{p) and the fermionic Vg{p) (up to a zero mode contribution). 
The zero mode as well as the doubly fermionic part of the free energy could then 
be found by rewriting Vg{p) as a total derivative in the fermionic potential. The 
purely bosonic part of the free energy can be calculated with the methods of ref. 
[8] . In principle the application of loop insertion operators to the free energy then 
yields arbitrary multi-superloop correlators. As we demonstrated for genus one, in 
practice these expressions become quite lengthy. 
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We then turned to the double scaling limit of the model in the moment descrip- 
tion. The m'th multi-critical point was identified and the scaling properties of the 
moments and basis functions were derived. The iterative scheme for the calculation 
of higher genus contributions could be optimized to produce only terms relevant in 
the double scaling limit. The general form of the free energy and the one-super loop 
correlators at genus g were stated. We presented explicit result up to genus two. 
Finally we commented on the identification of the supereigenvalue model at its m'th 
multi-critical point with 2d supergravity coupled to (4m, 2) minimal superconfor- 
mal matter. 

The analogy of structures in the Hermitian one matrix model and in the su- 
pereigenvalue model is rather impressive. This lends hope to find an answer to 
further interesting questions one should address, such as the supersymmetric gen- 
eralization of two- and multi-matrix models or the generalization of matrix models 
in external fields, which we shall be discussing in the next chapter. 



Chapter IV: 

The External Field Problem 



So far the partition functions of the Hcrmitian matrix and the supereigenvalue 
model acted as generating functionals at the same time due to the presence of 
arbitrary polynomial potentials. This is why solving for the free energy of these 
models really meant completely determining all possible correlators as well. As 
we know from field theory an alternative way to build generating functionals is to 
include external fields as sources coupled linearly to the fields of the theory. In the 
language of Hermitian matrix models this leads to partition functions depending on 
the eigenvalues of an external matrix. It is thus not surprising that interrelations 
between ordinary and external matrix models may be established. The key point 
here is to find a translation scheme from the coupling constants of the ordinary 
theory to the eigenvalues of the external matrices known as Miwa transformations. 
The study of matrix models in external fields has received much attention following 
the work of Kontsevich [44] whose external field model was shown to be equivalent 
to the Hermitian matrix model in the double scaling limit. No limiting procedure 
is necessary as the model directly produces only the scaling relevant terms. 

Moreover external field models are closely related to Hermitian two matrix mod- 
els, which in their scaling limit are capable of describing all the {p, q) minimal 
conformal field theories coupled to 2d gravity. Although the form of the two ma- 
trix model is obvious on the matrix level, it is not on the eigenvalue level as two 
arbitrary Hermitian matrices may not be simultaneously diagonalized by one U{n) 
transformation. 

In regard of these results it would be quite desirable to find the correct super- 
symmetric generalizations of external Hermitian matrix models. However, in that 
respect this chapter has the form of an outlook, as we have not been able to perform 
this program completely. Nevertheless some encouraging results may be reported. 

Due to simplicity we study the discrete case, i.e. away from the double scaling 
limit. After a review of the bosonic case we turn to the supersymmetrization of the 
observed structures. 

1. The Bosonic Case 

Let us consider a N x N Hermitian matrix model with a potential Vq{X) and an 
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external Hermitian matrix field L as source term 



Zn[L] = I VX exp[TriX - TrVblX)]. 

JnxN 



(4.1.1) 



By taking matrix derivatives of Zn [L] with respect to the source and defined by 
d 



{dL)ab 



dL 



ha 



{dL)abZN[L] - {Xab)ZN[L] 



(4.1.2) 



one may generate arbitrary correlators. In this sense the source L replaces the 
infinite set of coupling constants gk of the Hermitian matrix model of eq. (1.1.1). 

External field models may be solved in the planar limit and in their genus expan- 
sions through the method of Schwinger-Dyson equations [45, 46]. These equations 
are easily derived. Consider the shift of the matrix X in eq. (4.1.1) 



X ^ X + enX 



n+l 



n > -1, 



(4.1.3) 



with e„ infinitesimal. Then the invariance of the integral under a renaming of 
integration variables yields 



VXTt 



NxN 



L - Vo'(X) + ^ Tr [X'= ] 



k=0 



exp[ Tr i X - Tr Vo{X)]=Q. 



(4.1.4) 



in first order e„. 



□ 



See this by calculating the Jacobian of the transformation (4.1.3) to first order e„ 



Det 



dX' 



' Trln(l+.„ + ' ) Tre„a2£l!+i 

- = e oxT ' = e dxT = 1 . 



Now 



Tr ■ 



ax"+i _ a 
axT ~ axT 



.Tr 



dXT 
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IV: The External Field Problem 

and combining this with the variation of the potential and the source term yields eq. (4.1.4) 

Using derivatives with respect to the source eq. (4.1.4) may be rewritten as 



Tr 



L {dLr+' - V^idL) + ^(9^)'= Tr 1 Zn[L] = 0, (4.1.5) 



k=0 



and by pulling out all matrix derivatives we have 



Tiidi 



\n+l 



L-V^idL) Zn[L]=0 



(4.1.6) 



This system is in lact equivalent to a single matrix valued equation 

known as the Gross-Newman equation [45] . I — I 

Interestingly enough the external field model of eq. (4.1.1) with the potential 



(4.1.7) 



Vo{X) = -X^-n\nX 



(4.1.8) 



is closely related to the n x n Hcrmitian one matrix model with general polynomial 
potential 



[ PXexp[-Tr^gfeJ»Q 

Jnxn 



(4.1.9) 



The exact relation between the partition functions (4.1.1) and (4.1.9) was noted by 
Chekhov and Makeenko [47] and reads 



2n [ gk] = exp[- ^ Tr ] [L] , 



(4.1.10) 



provided L and gk are related by the Miwa-transformation 

for A: > 1 
2 ~ 

go = -Tr In L 



fffe = ^TrX-^--i4,2 



(4.1.11) 



in the limit N ^ oo which turns the g^-'s into independent variables. Note that 
in eq. (4.1.10) the N x N external field model is related to a n x n Hermitian one 
matrix model. The proof of eq. (4.1.10) is based on the fact that the Gross-Newman 
equation (4.1.7) for the potential (4.1.8) is equivalent to the Virasoro constraints of 
eq. (1.3.2) for the Hermitian one matrix model (4.1.9). 



4.1. The Bosonic Case 

Let us prove the relation (4.1.10). Taking the Gross-Newman equation (4.1.7) 
and applying one more matrix derivative to get rid of the singular term 
coming from the derivative of the logarithm we find 

dl-{n + l)-LdL)zN[L]=0, (4.1.12) 

which may be restated as 

[dl + LdL-n^ exp[-^TTL'']ZN[L] = 0. (4.1.13) 

Now insert the Miwa transformation into the matrix differential operator of this 
equation. By using the chain rule we have 



dL 



fc=0 



and 



(4.1.14) 



(4.1.15) 



fc=0 



Rewriting 



{dL 



oc -^ 



k,l 



k=0 a=l n=0 



tire double sum in the fi :st 

E E '^^ (Z°' ) 



n=0 

term of eq. (4.1.15) leads us to 

OO ^ 7n 



m— — 1 a— 1 m— k—0 

which is-already strongly reminiscent to the Virasoro constraints if one would ignore 



k=l 



fc=0 



{dl+LdL-n}^ J2 ;^{ETr(i^-')5,,„,,-a,„,.+Ey^.-.} 

(4.1.17) 



□ 



so by identifying k gk ^ Tr {L — 6k.2 in the above we recover the Virasoro gen- 
erator £jm 

OO ^ oc m 

{dl+Ldr.-n)= ;^{E'=5.-^.,„+.+E^^^^^."-> 



k=0 



k=0 



-{n + dg„), 



;4. 



1.18) 
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4.2. The Super-Miwa Transformations 



of eq. (1.3.3). Therefore the Vhasoro constraints obeyed by Zn[gk] are equivalent 
to the Gross-Newman equation for Zn [L] , provided that 



[gk] = -nZn[gk] 



(4.1.19) 



is vahd. But this is obviously true for an n x n Hermitian matrix model, as 
dgg Zn[gk] = ~(Ti'l) Zn[gk] seen from the definition in eq. (4.1.9). Hence we have 
proven the relation (4.1.10). Amazingly enough this identity does not depend on 
the dimension N of the external matrix field. 

Of course these identities may also be stated in the eigenvalue language. As the 
partition function of the external field model (4.1.1) only depends on the eigenvalues 
of the matrix source L, the integral over the angular variables may be performed 
with the help of the Itzykson-Zuber formula [48] to give the external eigenvalue 
model 



Z]sf[li 



(/n 



dx. 



T\i<j ("^i -^J ) 



N 

exp ^(-Vb( ) , (4.1.20) 



modulo an irrelevant 
values of X and Z,- 



multiplicative 
those of L. 



constant. The 
Integrating Zjv| 



denote 
with 



the eigen- 
the measure 



ni<j(^i — Ij)'^ exp[— J2i Viih)] over li yields the eigenvalue formulatlioil of the Her- 
mitian two matrix model. 

In terms of eigenvalues the Schwinger-Dyson equation (4.1.13) now takes the 
form 



1 



d_ 

dk 



exp[ 



\ZN[k 



d_ 

dk 



0, 



(4.1.21) 



k=l 



which may be derived by considering the variation 6li —eh in eq. (4.1.20) or by 
studying the eigenvalue dependence of the matrix derivatives Ol and 9£ . 

The identity (4.1.10) may also be proven directly on the level of th e ei genvalue 



models, by inserting the Miwa-transformations into eq. (4.1.9) and ufing integral 
identities for Hermite polynomials [49]. 

The equivalence of external field and Hermitian matrix models continues to hold 
true in the continuum limit. Here the Kontsevich model [44], which is built from 
a cubic Vo{L), is the matrix model describing the Hermitian matrix model in the 
double scaling limit, as shown in [50]. Let me just rr[en 
of generalized Kontsevich models a unification of all 
to 2d gravity could be achieved [51]. 



2. The Super-Miwa Transformations 

Our analysis of the supereigenvalue model in chapter HI revealed a striking num- 
ber of analogies to the Hermitian ma :ri: ; model. It is thus rather tempting to ask 
whether there exists an external field fo 'mulation of the supereigenvalue model as 
well, generalizing the results of the previous section. The partition function of 
such an external supereigenvalue model should depend on a set of even and odd 
external "supereigenvalues" h and 9i respectively. Similar to the Gross-Newman 
equation the Schwinger-Dyson equation of this model would then be represented 
as a differential operator in the external fields k and 6i annihilating the partition 
function. After performing a super-Miwa transformation of this differential opera- 



tor from the h and 6i to the coupling constants gk and ^, 



k+l/2 



one should recover 



the super- Virasoro generators, thus proving the equivalence of these two models. 

Unfortunately we have not succeeded in performing this program so far. However, 
we can report on some preliminary results. The correct differential operator in 
external fields could be identified, leading to the mjier- Virasoro constraints through 
a set of super-Miwa transformations. This operato: may even be formulated on the 
matrix level, thus giving hope to find a true supersymmetric external matrix model. 
Let us consider a pair ofr^ternal N x N matrices, the Grassmann even matrix L 
X A with Grassmann odd entfies. These shall be related 



accompanied by the matr 
to the coupling constants 5^ 
Miwa transformations 



md ^k+1/2 of the supereiger va 



1 



□ 



fffc-^TrZ-^- k>l 
go = -Tr Ini, 



efc+1/2 - -Tryli-^-^ fc>0 



ue mpd(fl via the super- 



(4.2.1) 



generalizing eq. (4.1.11). Note that at this stage we have not included delta function 
additions as in the bosonic case of eq. (4.1.11). The matrix derivatives 

' ° □ [3.2.2) 



d 



and 9yi 



may then be reexpressed as differential operators in the coupling constants by using 
the chain rule as 



dTvAL 



-k-l 



dAT 

k=0 feiO 



(4.2.3) 



and 



;ion, that in the framework 
q) minimal models coupled 



„ ^/l dTiL-'' „ dTiAL-''^^ „ 



fe=0 
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k=0 



fc+1 



J2(-L-'^-'d,,+J2L-^AL-'^-'+^d^,^,,, ). (4.2.4) 



As the super- Virasoro generators Gn+1/2 ^-nd iL„ are maximally quadratic in the 
dgi^ and S^^+i/a obvious that we should study quadratic matrix differentials in 
L and A. The first candidate is Da Ol which after some calculation may be found 
to be 



00 00 m 

Tr [dA 5i) = ^ Tr ( ) A: + 



fc=0 



fc=0 



fc=0 



Tr[dAdL+Adl] = 

00 oc 00 

m 
A:=0 

00 00 ^ (It 



Qm + k + l 



m— — 1 



fc=0 



fc=0 



00 ni~ I 



kd. 



m— — 1 



(4.2.5) 



E,, TO + 1 , . „ 
[k + ) Cfc+1/2 %+„+i/2 



k=0 
m — 1 



Note that already a number of factors present in the super- Virasoro generators of 
eq. (3.1.1) appear. The next candidate to consider is the expression Ad^. A rather 
tedious computation reveals that 



00 00 m 

Tr(Aai)= Tr(Ai-'"-2)[^fc^^a^^^^^ 

00 

^ k ^fc+1/2 



+ ^ Y fc^m-fc-l/2C'fc+l/2 } 
fc=0 

^ Tr(i-™-2)g^^^^^_,_2 Y Tr(Ai-™-2)L, 



(4.2.8) 



m— — 1 



fc=0 



fe=0 



(4.2.6) 



with the super-Virasoro generators Gm+1/2 and defined in eq. (3.1.1) ^. 

The relation (4.2.8) may also be stated on the level of "supereigenvalues" k and 
6i, despite the fact that we do not know how the 9i should be related to the matrix 
A. With the obvious supereigenvalue version of the super-Miwa transformations 



where all terms of order 3 in A cancel due to their fermionic structure. Additionally 
one proves the following identities 



00 00 
0= Y Tr(Ai— 2) [^(fc + TO + 1)^+1/2% 

m=-l k=0 



0= Y Tr(i— ^)[^e.+i/2a„„,,,, 

m=-l fc=0 



9k = lY^^-' k>l 

i 

go = "X!^^^*' 

i 

one has the foUowing identity 



u 



(4.2.9) 



d d d"^ y-^ 1 yd d a \ t ^ ^ 



E Tr(Ai- 



m— — 1 



k=0 



(4.2.7) 



E] ^ [ 5m+l/2 + 2 6*; /Cm ], 



(4.2.10) 



So by putting eq. (4.2.5), (4.2.6) and (4.2.7) together we see that the correct differ- 
ential operator yielding the super-Virasoro constraints is as simple as 



1 Note 


;ha 


; the factors of N'^ 


ik+1/2 N 




1/2- 
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as already noted by Huang and Zhang [52] in a different context. So we find ourselves 
in the situation where we know the matrix and eigenvalue based formulation of the 
Schwinger-Dyson equations of a model, but do not know the model they belong to. 
Including delta function additions in the super-Miwa transformations of eqs. (4.2.1) 
and (4.2.9) as they appear in the bosonic transformations of eq. (4.1.11) leads to 
terms linear in and ■ Unfortunately we have not succeeded in formulating the 
correct external supersymmetric model so far. 



3. Properties of External Eigenvalue Models 

In section 1 we saw how the Virasoro constraints emerged out of the Schwinger- 
Dyson equation associated with the shift SX = eX for the aseternal field model of 
eq. (4.1.1) with the special potential (4.1.8). There is yet mc ther way in which 
the Virasoro structure is hidden in the Schwinger-Dyson equSrions (4.1.6) of the 
external field models with arbitrary potentials as noted by Makeenko and Semenoff 



in ref. [46]. It turns out that this structure may be generalisec to an external 

supereigenvalue model. | 

Let us first study the bosonic case in the eigenvalue description. The external 
field model of eq. (4.1.20) reads 



Zn[1] 



N 



— 1 L J 



N 



(4.3.1) 



where A[Z] — Y[i<i{^i ^ the van der Monde determinant. The Schwinger- 



Dyson equations associated to the shift Xi - 
ten in the form 

L„A[/]Zjv[?] =0 



En for n > —1 may be writ- 



for 7{>|-1 



(4.3.2) 



with 



fc=0 ]=^i 



dl, 



Note that the generators L„ annihilate A[l] Zn[1]. One may easily construct gen- 
erators Cn which annihilate the partition function Zn [ I ] itself. For this purpose 
introduce the "long derivatives" 



(4.3.5) 



which commut' 
now read 

and 



i -with each other. The Virasoro constraints of eqs. (4.3.2) and (4.3.3) 
CnZN[l]=0 for n>~l (4.3.6) 



Cn = V/. ) (V,, )"+! + (V, k 

i 

1 " 

+9 EE(^'.)"(^'.r 



k=0 i=tj 



(4.3.7) 



Due to the commutativity of the V; . 's the generators of eq. (4.3.7) obey the Virasoro 
algebra (4.3.4) as well. 

For the supersymmetric case let us take the following ansatz for an external 
supereigenvalue model 



N 



Zn[1,1^]= j {J{dKde^^^^ex^^{-V{\) + Kk + e,^i,)\, (4.3.8) 

with the measure A"[A, ~ Y\i^j{\ ~ ^ ^i^j)"' ^-nd a undetermined. Note 
that the fermionic variables Oi enter only through the source term and the measure, 
can now study the Schwinger-Dyson equations of this model associated to the 
in integration variables 

<5A, = -e„ Ar+' and <5&, = e„ A,"+i (4.3.9) 

with n > — 1 and e„ an Grassmann odd parameter. The resulting correlator may 
be rewritten in the form 



(4.3.3) where 



Gr, 



The differential operators L„ obey the Virasoro algebra without central extension, 
as one verifies by direct computation ' — ' 



-1/2 ^ 



a,+i/2A"[Z,/i]Zjv[/,Ai] = for n>-l 



(4.3.10) 



[ Ln, ] — (n — m) i„+„ 



(4.3.4) 



□ 



d / d \ k / d \ "-k 



k=0 i^j 
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And in fact these generators obey the super- Virasoro algebra which one may directly 
check 



[ Ln,Lm ]^ (n-rnj Ln-i 



[ in, G„i+i/2 ] = Gn+m+1/2, (4.3.12) 

and where the generators L„ take the more complicated form 



iE(-i)(|-)VfEE(|-rfe 

n-1 



a 



In fact no constraints on the constant a in eq. (4.3.8) arise. In order to obtain super- 
Virasoro generators directly annihilating the partition function Z]\[[l,i2] introduce 
the "long derivatives" 

V„.A-.|,,H|-A"|i..l = |;-i:|^. («.14) 

Note that the Vi. commute and the V^. anticommute with each other. Replacing 
the derivatives di^ by V/. and the 9^. by V^^ in eqs. (4.3.11) and (4.3.13) then 
yields the super-Virasoro generators Gn+i/2 ^-nd which annihilate the partition 
function ZN[l,fi]. 

So we have seen that the "naive" ansatz of eq. (4.3.8) is capable of producing 
a generalization of the properties found in the external Hermitian matrix model. 
However all efforts to use this ansatz as the model generating the Schwinger-Dyson 
equation (4.2.10) associated to the super-Miwa transformations have failed. 

This concludes our short outlook on external supereigenvalue and supermatrix 
models. 
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